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C^ We review the present state of understanding of the Barkhausen effect in soft ferromagnetic materials. 
fjBarkhausen noise (BN) is generated by the discontinuous motion of magnetic domains as they interact 
q^with impurities and defects. BN is one of the very many examples of crackling noise, arising in a 

S'ariety of contexts with remarkably similar features, and occurring when a system responds in a 
erky manner to a smooth external forcing. Among all crackling system, we focus on BN, where a 
| * .c omplete and consistent picture emerges thanks to an exactly solvable model of avalanches dynamics, 
j^known as ABBM model, which ultimately describes the system in terms of a Langevin equation for 
+Jthe velocity of the avalanche front. Despite its simplicity the ABBM model is able to accurately 
^reproduce the phenomenology observed in the experiments on a large class of magnetic materials, 
■j_Jas long as universal properties are involved. To complete the picture and to understand the long- 
C^tanding discrepancy between the ABBM theory and the experiments — which otherwise agree 

axceptionally well — consisting in the puzzling asymmetric shape of the noise pulses, microscopic 
^ etails must be taken into account, namely the effects of eddy currents retardation. These effects 
i i y:an be incorporated in the model, and result, to a first order approximation, in a negative effective 
ijmass associated with the wall. The progress made in understanding BN is potentially relevant for 
pother crackling systems: on one hand, the ABBM model turns out to be a paradigmatic model for 
wthe universal behavior of avalanches dynamics; on the other hand, the microscopic explanation of 
Othe asymmetry in the noise pulses suggests that inertial effects may also be at the origin of pulses 
asymmetry observed in other crackling systems. 

7-H ■ 

> : 

^^L 1 Introduction 

t-H ■ 

Crackling noise 

^^Che term "crackling noise" [1,2] refers to the signal that some disordered sys- 
Q^tems produce as a response to an external driving field smoothly changing in 
Otime. Due to the presence of disorder, crackling signals are extremely irregular, 
^>clespite the steady increase of the external forcing. They are typically charac- 
I^jterized by a sequence of pulses of very different sizes and durations, separated 
;_jby quiescence intervals. Tiny events occur very frequently, while large ones are 
cdrare, with power laws probability distributions. 

In this review we focus on one particular case of noise, the so called 
Barkhausen noise (BN), which is emitted by ferromagnetic materials during 
the magnetization reversal process [3-5]. 
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Systems that "crackle" are found in many different situations, and, remark- 
ably, the corresponding signals often share some common characteristic fea- 
tures. Examples of crackling signals, besides BN, include the shear response of 
a granular media [6-8] , the acoustic emission during martensitic phase transi- 
tions [9], the bursts of dislocations activity in plastic deformation [10-14], the 
dynamic of superconductors [15-17] and superfluids [18,19], the fluctuations 
in the stock market [20-22], the dielectric polarization of ferroelectric materi- 
als [23], the acoustic emission in fractures [24-28], and the seismic activity in 
earthquakes [29,30]. 

Crackling noise signals are expected to encode information on the physical 
process that generates them. Understanding the statistical properties of these 
jerky emissions, is therefore a step towards the understanding of the micro- 
scopic dynamics taking place in the system that crackles. Moreover, the fact 
that very diverse systems behave in a remarkably similar manner, suggests 
that some general basic principle may exist in the underlying physics. If this is 
the case, then, understanding in detail one of them, may give insight to others, 
and propel other fields of research. 

BN can be seen as a case study, for which a satisfying picture is finally 
available. A main ingredient in this picture is a model of avalanche dynamics 
known as ABBM model, named after the authors of the original papers, Bruno 
Alessandro, Cinzia Beatrice, Giorgio Bertotti and Arianna Montorsi, whose 
companion seminal papers [31,32], one devoted to the developing the theory, 
the other to the comparison with experiments, have represented a huge step 
forward in the understanding of BN. The ABBM model, corresponds to a mean 
field description of an elastic magnetic wall moving in a disordered ferromagnet 
under the effect of an external driving field. The system is ultimately described 
by a simple Langevin equation for the velocity of the center of mass of the 
wall. 

Despite its simplicity, the ABBM model is able to reproduce with striking 
accuracy most of the phenomenology observed in BN experiments on a large 
class of magnetic materials. Moreover, thanks to a mapping onto a simple 
stochastic process, all the results concerning the noise statistics can be derived 
analytically, and have a clear and direct interpretation. In particular, this 
simple mean field description allows to explain the origin of the power law 
distributions of size and duration of pulses, both in the quasi-static limit, 
where the applied field variations are extremely slow, and at finite driving 
field rates. In this second case, the model is also able to predict how these 
distributions depend on the driving field rate. 

A more refined analysis of the experimental data against theoretical predic- 
tions includes the comparison of the shape of the avalanches. This is where 
the ABBM model unexpectedly fails, being unable to predict the characteris- 
tic leftward asymmetric form of BN pulses. The origin of this asymmetry lies 
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in the non-instantaneous response of the eddy fields to the domain wall dis- 
placement. To understand and evaluate the effect of such delay, one as to take 
into account the dynamical effect of eddy currents in the Maxwell equations 
for the eddy field. Eddy currents retardation gives rise to an "anti-inertial" 
effect, that can be accounted for, to the first order, by associating a negative 
effective mass to the wall. 

Once the effective mass is identified, the corresponding inertial term can be 
incorporated into the ABBM equation of motion. The modified ABBM model 
that comes out from this analysis gives a complete picture of the Barkhausen 
effect, consistent with the phenomenological observations: it is still success- 
ful in reproducing the distribution of sizes and durations of pulses, and also 
correctly reproduces the asymmetry in their shape. 

It is interesting to underline the role of universality versus non- universality 
in the theory of Barkhausen effect. The standard statistical mechanics ap- 
proach usually focuses on universal quantities. Definitely, universality is a key 
and extremely powerful concept at the basis of this approach, as it allows to 
predict the essential behavior of self-similar systems by means of very simple 
models. This consideration indeed applies to the case of Barkhausen noise, 
where, in analogy to critical phenomena, most of the statistical properties of 
the signals only depend on general properties of the system, while they are 
independent of the microscopic details: different magnetic samples respond to 
the forcing by an external field by producing events of magnetization reversal 
characterized by the same power law distributions, regardless of the specific 
microscopic structure of the material. 

However, some interesting features of BN turn out to be of microscopic 
origin. Although, consistently with universality, microscopic details would be- 
come negligible on extremely large scales, they have an unusually large effect 
in Barkhausen signals, an effect that is still significant at the experimental 
scale. On one hand this can be seen as an inconvenient and a limitation of 
the methods of statistical mechanics, since it means that not all the rele- 
vant aspects of the phenomenon can be captured by simple models. On the 
other hand, the identification of some macroscopic effect of the microscopic 
details, provides a mean to extract information on microscopic quantities by 
Barkhausen measurements. The asymmetry of the pulses gives indeed a mea- 
sure of this type. As an example, the skewness of the pulses measured as a 
function of the pulse duration shows a peak, that allows to identify a charac- 
teristic timescale for relaxation, which corresponds to the ratio between mass 
and damping constant. 

In this review we will deal with both these aspects of the theory, in the 
attempt to put together a number of results achieved in recent years. We will 
also try to underline which points are still unclear and would need further 
analysis. 
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1.2 Outline 

The paper is organized as follows. Section [2] reviews some general facts about 
Barkhausen noise. We start with a short historical note (|2.ip . Then we discuss 
how magnetic materials are classified, in order to clarify to which ones the 
following theory applies f|2.2j> . We shortly describe the typical set-up and the 
problems arising in BN experiments in subsection 12.31 BN phenomenology, as 
it comes out from experiments, is described in subsection 12.41 In subsection 
12.51 we discuss the two main theoretical approaches to BN that have been 
developed in recent times: one based on a microscopic description in terms of 
a random field Ising model, the other describing the magnetization process in 
terms of the dynamic of magnetic interfaces, that will be further developed in 
the next section. 

In section [3] we focus on the second approach: first we introduce the ABBM 
as a phenomenological model (13.11) and underline its main limitations (|3.2f) . 
To put the model in a more general framework we then briefly discuss the 
problem of the dynamic of a generic elastic interface in a random media (|3.3p , 
and then we specialize to the case of a magnetic interface (|3.4p . The purpose 
of this digression is to recover the ABBM model as a mean field for a model 
of interface dynamic f)3.5f) . therefore corroborating its validity. 

Section [4] is the core part of the paper: we show how the ABBM equation is 
mapped onto a Langevin equation describing a biased random walk in loga- 
rithmic potential with absorbing boundary conditions at the origin (|4.ip . This 
approach is amenable of analytical treatment and allows to work out exactly 
a number of properties, and to give a direct interpretation of the phenomenol- 
ogy observed in the experiments in terms of properties of the corresponding 
stochastic process ()4.2p . The power law distributions in avalanches' sizes and 
durations can be exactly calculated in terms of distributions of return times 
to the origin (|4.3p . The continuous dependence of the exponents in the power 
law distributions on the driving field rate can also be calculated, and turns 
out to be related to the marginality of the logarithmic perturbation to the free 
random walk (|4.4p . The existence of a threshold in the driving field rate can 
be derived from recurrence properties of a free random walk (|4.5p . and the 
existence of a cut-off in the power law distributions can be associated to the 
bias term (|4.6p . The average shape of the Barkhausen pulse can be derived in 
terms of the excursion of the process (|4.7p . and finally the power spectra can 
be derived from the process correlations (|4,8p . In subsection 14.91 we summarize 
the ABBM predictions and compare them with the phenomenological obser- 
vations. In the last subsection (I4,10|) we discuss the scaling of the cut-offs with 
the demagnetizing factor. This scaling is not entirely captured by the mean 
field approximation, although it is predicted by the full elastic interface model 
in d = 3. 
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Section 0] collects a number of results, most of which, but not all, are spread 
in literature, but to which has never been given a comprehensive exposition in 
the BN framework before. Some of these results have been indeed developed 
in an unrelated context, to deal with a completely different problem [33]. For 
this reason in this part we made an effort to give a self-contained exposition. 

The issue of the asymmetry of Barkhausen pulses (that the ABBM model is 
not able to reproduce) is treated in section [5j a characteristic leftward asym- 
metry in the average pulse shape is observed in experiments (|5.ip . which is 
due to the non-instantaneous response of the eddy currents field to the wall 
displacement ()5.2|) . The eddy current retardation gives rise to inertial effects, 
that are accounted for in a generalized ABBM model by associating a nega- 
tive effective mass to the domain wall (|5.3p . In subsection 15.41 we discuss how 
this non-universal effect is affected by the sample geometry, and, finally, in 
subsection 15.51 we compare theory predictions and experiments. 

In section [6] we discuss some other crackling systems that may be tackled 
with an approach similar to the one used for the Barkhausen effect: the case 
of the dynamic of a granular media under shear (|6. 1 1) . and the case of seismic 
activity during earthquakes (|6.2p . 

Section [7] and [8] are devoted to conclusions and acknowledgments respec- 
tively. 



2 The Barkhausen effect: history, experiments and theoretical 
approaches 

2.1 Short historical note 

Barkhausen noise is probably one of the first crackling signals ever recorded. 
Indeed the Barkhausen effect has been known for almost a century: its first 
observation dates back in 1919, when Heinrich Barkhausen noticed that "iron 
produces a noise when magnetized: as the magnetomotive force is smoothly 
varied [...] it generates irregular induction pulses in a coil wound around the 
sample that can be heard as a noise in a telephone" [3] . 

Barkhausen's first experiment consisted in observing the variations of the 
magnetization in a ferromagnetic sample subject to a slowly and smoothly 
changing magnetic field. A pickup coil winded around the sample was used 
to detect the variations of the magnetization as schematically shown in figure 
[TJ As the magnetization reverses, the variations of the magnetic flux induce a 
voltage in the coil, that indirectly measures the changes in the magnetization 
of the sample, and that may eventually be translated in real noise through an 
amplifier. The signal recorded in correspondence to jumps in the magnetiza- 
tion, appears to be very irregular, no matter how smooth is the variation of 
the external field. 



February 18, 2009 



15:51 



Advances in Physics 



adpl.21 



6 



Francesca Colaiori 



pick-up coil 



\ 



ferromagnetic sample 





induced voltage 



Figure 1. Sketch of experimental set-up used in the first BN experiment. The bar is the 
ferromagnetic sample, which is inserted in the pick-up coil. As the external field is changed, an 
electromotive force is induced in the coil. The voltage can eventually be transformed in real noise 

through an amplifier. 

Nowadays, magneto-optical methods are also available, that allows a direct 
observations of moving domains, however, the most reliable experiments in 
terms of signal statistics are still based on inductive measurements. In order 
to be capable to record clean signals with reduced background noise and high 
statistic, modern experimental settings make use of sophisticated equipments, 
however, the basic set-up for inductive BN experiments is conceptually still 
the same as the one shown in figure [TJ 

Barkhausen wondered about what could be the source of the signal emitted 
during the magnetization of the iron sample, and he mistakenly concluded that 
the sequence of pulses was generated by a corresponding sequence of sudden 
and complete reversals of entire magnetic domains. At that time, his experi- 
ment was indeed taken as the first indirect observation of magnetic domains, 
whose existence had been postulated a few years before (1907) by Weiss [34]. 
Almost 20 years later (1938) Elmore was able to perform an experiment on a 
cobalt crystal where the motion of domain boundaries was directly observed 
and measured for the first time [35]. However Elmore did not recognize this 
motion as the source of BN. Only more then 10 years later (1949), after an- 
other experiment by Williams and Shockley on a FeSi crystal [36], the origin 
of BN was finally ascribed to the correct cause, namely to the irregular fluc- 
tuations of magnetic domain boundaries, rather then to the sudden inversion 
of domains. This fact was definitely clarified in the same year by Kittle in his 
fundamental review on the physical theory of ferromagnetic domains [37]. A 
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nice exposition of the first steps in the understanding and misunderstandings 
of Barkhausen effect is given in a recent review by G. Durin and S. Zapperi [5], 
which also includes in appendix a translation from German to English of the 
original paper by H. Barkhausen. 

Since its first observation, a lot of work has been done in order to "decode" 
the Barkhausen signal, both on the experimental side and on the theoretical 
one. For a long time, theoretical studies have been based on a phenomeno- 
logical approach, which described the signal as a superposition of random 
elementary jumps. Then, it became clear that BN could represent a powerful 
tool to investigate the magnetization process and the hysteretic properties on 
a microscopic scale, and much effort has been devoted in developing physically 
grounded models that could allow to relate the phenomenology to its micro- 
scopic origin. The theoretical studies have followed two kind of approaches: 
one based on a microscopic description in terms of spins in a random magnetic 
field (see section [275]) ; the other one, which is the one on which we will focus in 
this review, describes the system in terms of a fluctuating magnetic interface 
that moves under the action of the external field. 

The research on BN has been in part motivated, especially at the begin- 
ning, by the applications of Barkhausen effect to material testing. Barkhausen 
emissions are indeed commonly used to check in a non-destructive way the 
integrity of magnetic samples: the signal intensity is sensitive to the changes 
in material microstructure, and to the presence of residual stresses (in mag- 
netostrictive positive materials compressive stresses will decrease the intensity 
of Barkhausen noise while tensile stresses will increase it). These properties 
make BN an efficient tool for the detection of micro-imperfections and for the 
evaluation and mapping of the local distribution of residual stresses [38-45] . 

Most of the recent research on the statistical properties of BN is however 
theoretically oriented, and primarily aimed to gain understanding on the hys- 
teretic properties of ferromagnetic materials, and to investigate the magneti- 
zation reversal process on a microscopic scale. 

Recently, further interest has been raised by the fact that many other driven 
dissipative systems have been found to respond to a smooth forcing in a dis- 
ordered manner, with very similar and reproducible statistical features. To 
this kind of response has been given the general denomination of "crackling 
noise" . Clearly, the irregular response is due to the presence of disorder in the 
system, and the challenge is to get information on the specific system under 
study by decoding the noise that it generates. As the recent encouraging ad- 
vances in the understanding of BN compose a rather satisfying picture of the 
phenomenon, BN emerges as an attractive case study to understand crackling 
noise in general. 
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Figure 2. Sketch of a typical hysteresis loop for a hard magnet (a) and a soft magnet (b). 

2.2 Soft magnets versus hard magnets 

The understanding of hysteretic properties of magnets has an interesting over- 
lap with material science, and the theory that we describe has been built al- 
ways keeping close contact with the experimental side. Therefore, although 
the aim of this review is mainly theoretical, we believe that, in order to put in 
a correct framework the theory, understand its limitations, and to keep close 
the connection with the experimental counterpart, it is worth to discuss briefly 
how magnetic materials are classified, and to clarify to which ones this theory 
applies, and to which ones it does not. 

Magnetic materials are roughly classified in soft and hard. Soft magnets 
are those characterized by very tiny hysteresis loop, reflecting the fact that 
they are easily magnetized and demagnetized by applying relatively small 
fields. Hard magnets are the exact opposite: they have very wide and squared 
hysteresis loops, and they need a big effort in terms of applied field in order to 
be magnetized and demagnetized. The two classes of magnets are quite well 
separated, and there is not much in between. A good quantitative measure of 
"magnetic hardness", is given by the coercive field, which is the field needed 
to bring back to zero the magnetization of a saturated sample. This quantity 
differs, between soft and hard magnets, by several orders of magnitudes (see 
figure ED [46]. 

The markedly different hysteretic behavior observed in the two classes of 
materials is due to the fact that the magnetization reversal proceeds by differ- 
ent mechanisms in the two cases. In soft magnets it is dominated, except than 
close to saturation, by the nucleation and displacement of magnetic walls. 
This is a process with low energetic cost, that translates into a small loop 
area, which measures the dissipated energy. On the other side, in hard mag- 
nets the wall nucleation and/or displacement is strongly impeded, or totally 
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suppressed, and magnetization reversal by coherent spin rotation enters into 
play. 

Hard magnets have highly anisotropic magnetic properties. They are usually 
composed from alloys containing rare earth elements (typical ones are Co^Sm 
and Ne2Fe\4 : B). In single domain particle magnets the nucleation of domains 
is prevented by building samples made out of very small magnetic grains, 
whose size is comparable with the domain wall thickness, and each one with 
a random easy magnetization direction. Due to the high anisotropy, the mag- 
netization can only reverse by coherent spin rotation when the applied field is 
large enough to force the spins out of the easy magnetization direction. Due 
to the high coercive field, hard magnets cannot be accidentally demagnetized 
when subjected to small fields. Also, they have a remanent magnetization very 
close to the saturation value, meaning that they are able to retain their mag- 
netization when the applied field is removed. These properties are extremely 
desirable in the design of magnetic storage devices. 

Examples of soft materials are Iron-Silicon alloys, crystalline Iron-Nickel 
alloys (permalloys), and amorphous metallic alloys (typical is the FeCoB). 
Applications of soft magnetic materials exploit the large flux changes obtained 
with a small change in the applied field, and the small energy loss. Soft mag- 
nets are used for example in transformers, motors, and inductors, and also as 
field sensors in magnetic recording devices. The magnetization reversal in these 
materials is dominated by the low-cost process of magnetic domain displace- 
ment. Domain walls, and the domain walls movement can be directly observed 
in experiments based on the magneto-optical Kerr effect (MOKE). Reflection 
of a beam of linearly polarized light from the surface of a magnetized sample 
causes the polarization to become elliptical, with the principal axis rotated 
with respect to the incoming light, depending on the magnetization direction. 
An analyzing filter in the reflected beam is then used to generate magnetic 
domain contrast. 

The theory described in this review is based on a mesoscopic description 
in terms of dynamics of domain walls, therefore it only applies to soft mag- 
nets. The magnetization reversal process in hard magnets is indeed worth of 
attention, however, in this case, a microscopic description in terms of spin 
models would be more appropriate. As it will be discussed in subsection 12.41 
soft materials have different statistical properties depending on the range of 
the interactions ruling the dynamics. The mean field theory which is the focus 
of this review applies to those with long range interactions. 

2.3 Barkhausen noise experiments 

In inductive experiments BN is detected by a pick-up coil wound around the 
ferromagnetic sample during the magnetization of the material, which takes 
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place under the action of a varying external field H(t). As the magnetiza- 
tion reverses in the sample, the variations of the magnetic flux <E>, induce an 
electromotive force V in the pick-up coil, which has a contribution from the 
applied field across the pick-up coil, and another one given by the magnetic 
field inside the material: 

V = -iV$ = -Nfi (A C H + AM) , (1) 

where /io is the vacuum magnetic permeability, N is the number of coil turns 
A c is the coil cross section, A is the sample cross section, and M the magne- 
tization. 

In high susceptibility materials the contribution coming from the first term 
in equation ([I]) can be neglected, since, as long as the susceptibility % = 
(/Va*0 — 1) is much larger then 1 (where [i is the magnetic permeability of 
the ferromagnetic material), the contribution from the magnetization rate M 
is much larger then the one due to the variation of the applied field H. To 
neglect the first term in equation ([T|) , it is also necessary that the sample cross 
section A is not too small with respect to the pick-up coil section A c . 

In low susceptibility materials, or in cases where the condition A — A c is 
not easy to satisfy, as for example in the case of thin films, it is possible to 
compensate the induced flux in air by using another pick-up coil with the 
same number of turns and cross section, but wound in the opposite direction. 

With the eventual addition of the compensation coil, it can therefore al- 
ways be assumed that the variation of the magnetic flux <3? is approximately 
proportional to the rate of change of the magnetization: 

V ~ -Nfi AM . (2) 

In the simplest experimental conditions, where there is just one domain wall 
separating two regions of opposite magnetization, the rate of change of M is 
twice the magnetization of the region spanned per unit time. The volume of 
this region is given by vdL in the geometry and notations of figure [3j and 
where v is the average wall displacement velocity. This volume, times the 
magnetization per unit volume M S /AL, gives the magnetization rate 

M = 2M s vd/A , (3) 

where M s is the saturation magnetization. 

Therefore, the induced voltage V measured by the pick-up coil is just pro- 
portional to the wall velocity v: 



V = -2Nfi dM s v . 



(4) 
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Figure 3. Sample geometry for a typical inductive BN experiment. 



Note that the sample thickness d in the direction perpendicular to the wall 
movement and to the coil axis is the only geometrical parameter that enters 
in the intensity of the BN signal. 

The number of turns in the pick-up coil amplifies the signal by a factor N, 
thus one could be tempted to use wide detection coils. However, we stress that 
this is not a good choice: the magnetostatic field produces a counterfield which 
has a crucial role in the domain wall dynamics. This counterfield depends on 
sample geometry and domain structure, and it is homogeneous only in special 
geometries. However it is convenient to assume it to be constant also in more 
general cases. This is a good approximation only as long as one limits the 
region of the pick-up coil. 

In the more complex case where several domain walls are present, the signal 
encodes the complex collective effect of many interacting domains. In this case 
the quantity v in equation (J3|) approximately measure the velocity of the active 
walls. 

Due to the presence of various types of disorder, like non-magnetic inclu- 
sions and dislocations, the wall movement is discontinuous, since the disorder 
locations act as pinning points. As the external field is smoothly increased, 
the magnetization changes in steps, in correspondence to jumps of the mag- 
netic interface, or, in microscopic terms, in correspondence to avalanches of 
spin flips. The hysteresis loops is also discontinuous, and the signal, which is 
proportional to the derivative of the magnetization, plotted versus time, looks 
like a disordered series of pulses (see figure E|) . 

Note that, in identifying the Barkhausen signal with the derivative of the 
hysteresis curve, we are assuming that the applied field H(t) is increased lin- 
early in time: indeed, in BN measurements, often — but not always — a driving 
field with a triangular time profile is applied, instead of one with the sinusoidal 
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Figure 4. Sketch of a typical hysteresis loop for a soft magnet. The loop is not smooth, rather the 
magnetization increases in jumps as shown in the zoom. The inset on the bottom right side is the 
derivative of the magnetization, proportional to the typical Barkhausen signal, detected from the 
induced voltage. The voltage directly measures the velocity of the domain wall: pulses correspond 
to jumps of the magnetic interface that gets repeatedly stuck by the various kind of impurities and 

imperfection, acting as pinning centers. 

shape commonly used in hysteresis experiments. In this case, the variations 
of magnetization with respect to the applied field are directly translated in 
derivatives with respect to time. 

There exists a vast literature reporting results and data about BN experi- 
ments in soft materials. However, not all the data collected are reliable when 
one is interested in statistically analyzing and modeling the self-similar prop- 
erties of the noise, as the details of the experimental procedure may have a 
serious impact on the quality of the signal. Unfortunately, the need to establish 
a standard in the experimental set-up that could allow a rigorous comparison 
of the data, has been recognized as an important issue only in recent times. 

The most relevant factor to take into account in order to ensure a reliable 
noise signal is to ascertain the stationarity during the dynamics. This question 
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Figure 5. Hysteresis loop of a FegsBis amorphous alloy under moderate tensile stress (lOMPa). 
The signal measured along the whole half cycle (from A to A') is highly non— stationary. The 
interval BB' represent the interval where the permeability can be assumed to be constant, and 
where the Barkhausen signal is stationary. (Reprinted with permission from [48]). [48] Copyright 

2006 by Springer. 

is crucial, since non-stationarity strongly biases the scaling of the avalanches 
distributions. However, it is not always easy to check whether the system is in 
a stationary state or not. 

In the context of BN, the statistical properties of the signal strongly depend 
on the region of the hysteresis loop where the measure is recorded. This ques- 
tion was first raised by Bertotti and Fiorillo [47], who suggested to restrict 
the measures to the part of the loop corresponding to constant average per- 
meability, where the condition of stationarity is guaranteed. This important 
point has however been disregarded in many experimental studies. 

When collecting BN data, it is indeed crucial to take into account that the 
Barkhausen signal is stationary only in the linear region around the hysteresis 
loop, where the intensity of the applied field is close to the coercive field value 
(see figure [5]). This is precisely the region where the domain wall displacement 
is the dominating — if not the unique — mechanism involved in the magneti- 
zation reversal process. Getting closer to saturation, other mechanisms enter 
into play, as the spins, which during the previous phase of the magnetization 
process stay parallel to the easy magnetization axis, start to coherently rotate 
in order to align to the external field (see figure [6]). 

Experimental studies that do not take this point into account, and aver- 
age Barkhausen signals over the whole hysteresis loop, will mix the effect of 
the different magnetization mechanisms and of different dynamical conditions, 
making very difficult to attribute to the corresponding data a clear physical 
interpretation. 

The role of stationarity has been recently further analyzed in great detail 
by Durin and Zapperi [48]. In their work they show both with simulations and 
experiments that in non-stationary conditions the distributions of avalanches, 
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Figure 6. Typical magnetization loop for a soft magnet. On the right hand side the magnetization 

domain structures corresponding to five different regions of the hysteresis loop are sketched, 
assuming that only a single domain wall is present. The domain structure in (a) corresponds to a 
zero magnetization configuration. The domain wall divides the sample in two regions of opposite 
magnetization, oriented along the material easy magnetization directions. The two small triangular 
domains are the so called closure domains. Their presence is neglected in the following treatment. 

Increasing the applied field, the domain wall moves as in (b) in order to increase the size of the 
domain magnetized closer to the direction of the applied field. In this region the hysteresis loop is 
linear. This motion proceeds as the external field is increased, until the wall has spanned the whole 
sample as in configuration (c). At this point the sample is all magnetized in the same direction, 
which however does not coincides in general with the direction of the external applied field. From 
now on, further increasing the field, the magnetization proceeds by coherent rotation of the spins 
(d), until the saturation magnetization is reached (e), where all the spins are all aligned parallel to 

the direction of the applied field. 



being integrated over different values of the control parameters, are charac- 
terized by larger effective exponents. This might explain some experimental 
results which give unusually large values for the exponents in the power law 
distributions of sizes and durations of avalanches [49]. Evaluating the effects 
of non-stationary conditions on the statistical behavior of a system is an issue 
that arises in many driven non-equilibrium situations, and is analyzed in a 
more general context in [50]. 

A detailed analysis of the data from many Barkhausen experiments reported 
in literature is given by Durin and Zapperi in their review [5], where the 
reliability of each experiment considered is also discussed. 
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Figure 7. (a) Amplitude probability distribution Po(4) (which corresponds to P(V) in our 
notations) at three different values of the driving field rate. <j> is the variation of magnetic flux 
assumed to be proportional to the wall velocity, c is proportional to c. The experimental data on a 
polycrystalline 3% FeSi sample are fitted with equation [5] (b) Same as (a), in semilogarithmic 
scale. (Reprinted with permission from [32]). [32] Copyright 1990 by the American Institute of 

Physics. 

2.4 Barkhausen noise phenomenology 

BN has some remarkable statistical properties, that are nowadays measured 
with great precision in experiments. We briefly report the main experimental 
results here. 

2.4.1 Signal amplitude distribution. The amplitude V of the BN signal, i. 
e. the instantaneous value of the voltage, is found to follow a distribution well 
fitted by 

P(V) oc y-(!- c ) exp (- V/Vb) , (5) 

where c is proportional to the driving field rate, and Vo is some characteristic 
cut-off value (see figure [7]). 

The amplitude distribution of the signal is an average quantity over 
avalanches. To go more into detail, and investigate the statistical properties 
of the single avalanches, it is necessary to define precisely where a jump starts 
and ends. This is commonly done by fixing a threshold value, which cannot be 
strictly zero due to the unavoidable presence of background noise (see figures 
[U [9]). That of signal thresholding is however a very complicated issue, which 
arises in many cases when one has to analyze temporal series, and that would 
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Figure 8. Size of an avalanche. 



duration 




Figure 9. Duration of an avalanche. 

need further attention. The choice of the threshold is crucial, since it defines 
the objects whose statistical properties one is looking at, and the absence of a 
clear rule makes the identification of single pulses rather ambiguous. On the 
other hand, signal thresholding is inescapable if one wants to analyze directly 
experimental time series, it would therefore be desirable to be able at least 
to control its effects on the measured quantities. An alternative approach to 
the statistical analysis of the signals, which bypasses completely ambiguous 
thresholding methods, is to focus on the spectral properties, which indicate 
temporal correlations in the system, and then try to infer from there the statis- 
tical properties of the avalanches. This is indeed possible in some cases thanks 
to the fact that the power spectra exponents are related to the exponents in the 
size and time distribution of pulses as it was shown by Kuntz and Sethna [51] 
for the random field Ising model and some of its variants, and by Laurson, 
Alava, and Zapperi [52] for the case of sand piles. 

In the case of BN, the experimental data have been shown to be relatively in- 
sensitive to the exact choice of the threshold [53] . Background noise limits the 
detection of small avalanches. Choosing the threshold below its typical ampli- 
tude would result in over-estimating very small avalanches, but it is expected 
not to affect sensibly the rest of their distribution. The characteristic signal 
amplitude Vq in equation ([5]) gives a good reference value for the threshold, 
which is usually taken between 5 and 15% of Vq. 
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2.4.2 Distribution of sizes and durations of Barkhausen pulses. Once the 
pulses are well defined one can look at the statistics of their sizes and durations. 
By drawing an histogram it appears that they are also both characterized by 
power law distributions, limited by an exponential cut-off: the probability of 
having a pulse of a given size has the form 



where F$ is some scaling function that accounts for the cut-off, and So is some 
characteristic avalanche size. 

The distribution of durations has a similar behavior: the probability of hav- 
ing a pulse of duration T is given by 



where again To is some characteristic time, and Ft accounts for the exponential 
cut-off. 



2.4.3 Universality classes. By repeating BN experiments on different sam- 
ples of different materials and in different external conditions, it has been 
observed that the exponents characterizing the power law distributions are 
universal, and tend to cluster around two sets of values (see figure [10]) [54]. In 
soft magnets, this allows to identify two distinct universality classes, which are 
named long range and short range, according to the kind of elastic interaction 
which dominates the behavior of the magnetic walls. 

The two classes are characterized by exponents that, in the zero frequency 
limit of the applied field, take the values r ~ 1.50, a ~ 2.0 for materials in 
the long range class, and r ~ 1.27, a ~ 1.5 for those in the short range class, 
respectively for the avalanches' sizes and durations distributions (see figure 



Typical materials belonging to the long range class are poly crystalline 
FeSi with high Si content, and partially crystallized Fe§^Co2\B\§ amor- 
phous alloys. Short range behavior is observed for instance on Perminvar 
(_Fe3o-/Vi45Co25) and amorphous alloys with composition Fe x Co^- x B\^ under 
tensile stress. 

We will restrict here to materials belonging to the long range class, since for 
these systems the mean field description given by the ABBM model, which is 
the focus of this review, turns out to be exact for universal properties at long 
length and time scales (see section 13.51) . 



P(S) = S- T F S (S/S ) 



(6) 




(7) 



nop. 
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Figure 10. The data refer to Barkhausen measurements in different materials at the lowest 
available driving frequency, (a) Distribution of size of avalanches for different long range (filled 
symbols) and short range (empty symbols) materials. A solid line corresponding to slope 1.5 (long 
range) and a dashed line corresponding to slope 1.27 (short range) are shown for comparison, (b) 
Distribution of duration of avalanches for different long range (filled symbols) and short range 
(filled symbols) materials. A solid line corresponding to slope 2.0 (long range) and a dashed line 
corresponding to slope 1.5 (short range) are shown for comparison. (Reprinted with permission 
from [54]). [54] Copyright 2000 by the American Physical Society. 



2.4.4 Distributions of pulses at finite driving rates. The values r ~ 1.5 
and a ~ 2.0 for the long range elasticity class, refer to experiments performed 
in the quasi-static limit, namely, such that the applied field is increased very 
slowly with respect to the typical time scale of the system, in such a way that 
the system has time to rearrange after every infinitesimal change of the field 
before the field is increased again. 
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Figure 11. Cumulated distribution of sizes P(S) = p(s)ds for a polycrystalline FeSi material 

at different applied field frequencies. The solid lines are fits with equation [6] and with 
r — 1 = (1 — c)/2. The exponent r varies linearly with frequency between 1.5 and 1.15. (Reprinted 
with permission from [5]). [5] Copyright 2005 by Academic Press. 

In experiments where the frequency of the field is finite but still small, one 
can observe a linear dependence of the exponents a and r on the driving field 
rate c, well fitted by the expressions a = 2 — c, and r = (3 — c)/2 (see figure 
[TT| where the cumulated distribution of sizes of BN pulse for a polycrystalline 
FeSi sample is shown for seven values of the field sweep rate). No similar 
dependence has been observed in materials belonging to the short range class, 
where the exponents keep their quasistatic limit values, independently of the 
driving field rate. 

2.4.5 Threshold in the field rate, transition from intermittent to contin- 
uous regime. The variation of the exponents a and r with the driving field 
rate c reflects the fact that, increasing the sweep rate, large events become 
more frequent. Indeed, even without a statistical analysis, it is clear by eye 
inspection that the Barkhausen signals become more and more pronounced as 
the frequency of the applied field is increased (see figure PT2j) . 

Further increasing the driving field rate, it is observed that, above some 
well defined threshold, the typical intermittent behavior is lost, single pulses 
disappear, and the magnetization reverses in a unique avalanche. This is a sort 
of meta-phase transition, which separates two completely different regimes of 
the dynamics. It would be interesting to see how the system approaches the 
critical driving rate, by monitoring for example the quiescence times separating 
the avalanches as c approaches the threshold value. This, to our knowledge, 
has never been done. A similar situation arises for example in sheared granular 
matter, where a threshold driving separates the stick-slip form the sliding 
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Figure 12. Barkhausen signal recorded at increasing sweep rates. (Reprinted with permission 
from [5]). [5] Copyright 2005 by Academic Press. 



regime. This effect is due to avalanches overlapping aver time: since at finite 
driving rates the external field increases during the avalanche, it can keep 
active an avalanche that would have stopped otherwise (see subsection W. 
and the paper by White and Dahmen [55]). 



2.4.6 Cut off in the distribution of sizes and durations. The power law 
distributions of sizes and durations of avalanches are limited by an exponential 
cut-off, which turns out to be related to the demagnetizing field (see subsection 
I4.6|) . which in its turn depends on the sample geometry. Experiments repeated 
on samples with different geometries, indeed confirm that the position of the 
exponential cut-offs So and To in the distributions §§§ varies with the shape of 
the sample. However, a precise experimental evaluation of the cut-off depen- 
dence on the parameter k measuring the strength of the demagnetizing field 
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Figure 13. Top graph: experimental estimations of the cut— off So of Barkhausen pulse size 
distribution as a function of the demagnetizing factor k in a FeSi 6.5 wt% alloy (empty symbols), 

and an amorphous alloy Fe2iCo§^B\5 under stress (filled symbols). The empty symbols 
correspond to the long range class considered here. Bottom graph: same as the top for the cut-off 
To of the Barkhausen pulse duration distribution. (Reprinted with permission from [54]). [54] 
Copyright 2000 by the American Physical Society. 



is very difficult, also due to the further dependence of the distribution on the 
driving rate. 

Earlier studies [56] report the scaling So ~ fe _1 , and To ~ k~~ 1 / 2 . The most 
reliable experiments are probably those reported in [54,57], done by taking 
measurements on the same ribbon repeatedly cut along one direction, in order 
to modify the effect of the demagnetizing field, while stresses and internal 
disorder are kept constant. For a FeSi sample, the scalings reported are So ~ 
k-°- 57 , and T ~ A;- 30 (see figure Q3). 



2.4.7 Average pulse shape. Averaging avalanches of the same duration T, 
one obtains, for every given T a characteristic pulse shape. 
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t/T 

Figure 14. Attempt to collapse average pulses of different durations onto the same function. The 
data are from a BN experiment on a polycrystalline FeSi sample. Pulses are normalized with 
N oc Jg 1 d{t/T)(v{t | T)). The curves are averages over several pulses for each of eight different 
values of T, ranging from 0, 16 to 2, 03 ms. (Reprinted with permission from [5]). [5] Copyright 

2005 by Academic Press. 



(V(t I T)> (8) 

In analogy to conventional critical phenomena one might expect that average 
pulses of different durations, properly rescaled by some power of their duration, 
will depend on t only through the scaling variable t/T, and therefore collapse 
onto the same universal function. This quantity was recently proposed in [1] to 
test theory against experiments, and indeed it provides a much stringent tool 
than the simple comparison of scalar quantities, such as critical exponents. 

As shown in figure [14"1 which reports data from a BN experiment on a 
polycrystalline FeSi sample, the rescaling of experimental pulses can approx- 
imately be performed. As we will see, the shape of the scaling function is 
influenced by some non-universal effects that however disappear, as expected, 
on very long time scales. 

2.4.8 Power spectra. In most of earlier literature, experimental papers fo- 
cus on the detection of the power spectra of the noise. The most accurate 
report of power spectra are reported in the seminal paper by Alessandro et 
al. [32] in experiments on a 3% FeSi sample. For experiments on single crystals 
and other polycrystalline FeSi materials see also [58]. 
For materials in the long range universality class the power spectra is found 
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Figure 15. Log— log plot of the power spectra divided by the driving field rate, versus frequency 
measured in a polycrystalline 3% FeSi sample. The three solid lines are the experimental results 

for three different driving field rates, the dotted lines are predictions from the ABBM model. 
(Reprinted with permission from [32]). [32] Copyright 1990 by the American Institute of Physics. 

to decay at large frequencies uj as u>~ 2 (see figure [To]) , although in same cases 
a deviation from this scaling is observed at intermediate frequencies. At small 
frequencies the spectrum is characterized by the presence of a maximum. The 
position of the peak depends on the driving field rate c: it occurs at a frequency 
com roughly proportional to c 1 / 2 , and with an amplitude Fm that scales as c 
in the intermittent regime (slow driving), while it is roughly constant in the 
continuous regime (fast driving) (see figure [16] where Fm/c is plotted versus c). 

The behavior of the power spectra reported in [32] and shown in figure [15] 
is indeed common to many ferromagnetic materials. Similar spectra are found 
for example in NiFe alloys with vanishing anisotropics [59], grain oriented 
FeSi alloys [47], FeSi single crystals [60]. 

At frequencies below the peak the power spectrum has been claimed to 
behave like some power <J^ , with approximately if} ~ 0.6. However, these 
claims are based on fits of data that are spread on less than one decade, and 
therefore are not conclusive. 

Although the large u tail and the existence of the peak are reproduced by the 
theory, the correct interpretation of all the details in the complex behavior of 
the frequency spectrum still represents a challenging problem. In principle the 
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Figure 16. (a) Amplitude of the maximum Fm in the frequency spectrum divided by driving field 

rate versus driving field rate in log— log scale. Different data-sets refer to different values of the 
permeability. The bold line separates the two regimes intermittent and continuous, c is equivalent 
to c in our notations, (b) Same as (a), for the frequency u>m corresponding to the maximum of the 
power spectrum. (Reprinted with permission from [32]). [32] Copyright 1990 by the American 

Institute of Physics. 

contributions from the correlations among different avalanches, and that from 
the correlations within avalanches are mixed together in the power spectrum, 
which therefore is not easily amenable of physical interpretation. 

Some authors have also considered higher order spectra to further charac- 
terize the BN signal [61,62]. In principle this analysis could indicate whether 
events on a certain scale would systematically precede events on another scale. 
While no precursor events have been found, time asymmetries of the signal 
emerge from the analysis of the imaginary part of the moment of order 1.5 
consistently with that observed in the skewed pulse shape (see below). 

2.4.9 Summary of experimental observations. It is useful at this point to 
summarize the main phenomenological observations collected from different 
BN experiments, that we will want to compare to the theoretical predictions 
in the next sections. 

For soft materials with long range interactions: 

(i) The distribution of the intensity of the BN signal follows a power law decay 
with an exponent 1 — c where c is proportional to the driving field rate. 

(ii) The sizes and durations of avalanches are power law distributed. In exper- 
iments performed in the quasi-static limit the exponents are a ~ 1.5 and 
r ~ 2.0, respectively. 

(iii) The distributions of sizes and durations of avalanches flatten as the sweep 
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rate is increased. The corresponding exponents vary linearly with the field 
rate c, and their values are well fitted by the relations a = 2 — c, and 
T = (3-c)/2. 

(iv) The previous statement hold as long as the sweep rate is not too high: there 
exist a well defined threshold in the driving field rate c above which the 
Barkhausen pulses become indistinguishable and magnetization reverses in 
a unique avalanche. 

(v) The exponential cut-offs So, and To in the distributions of sizes and du- 
rations of avalanches are related to the sample geometry through the de- 
magnetizing field, and scale as powers of k. 

(vi) Average pulses properly rescaled by their duration can approximately be 
collapsed onto the same function. 

(vii) The power spectrum decays as uj~ 2 at large frequencies. At small frequen- 
cies is characterized by a peak that occurs at u>m — c 1 / 2 , with an amplitude 
that in the slow driving regime scales roughly as Fm — c. 



2.5 Modern theoretical approaches to Barkhausen noise 

Recent theoretical studies on BN have moved from a pure phenomenological 
description, where BN is considered as a superposition of elementary events, 
without any connection to the microscopic dynamic, to more physical and de- 
tailed descriptions which take into account the relevant interactions governing 
the dynamics of the magnetization reversal process. 

The theory has then developed on two different but related avenues: one 
follows a microscopic description and models the magnetic system through 
interacting spins. In this framework, the criticality observed in BN experiment 
is explained in terms of the proximity of a disorder driven non-equilibrium 
phase transition, which is observed in the zero temperature random field Ising 
model (RFIM). The other approach describes the magnetization process in 
mesoscopic terms, through the dynamics of the domain walls. In this case, the 
criticality is related to the depinning transition occurring when the applied 
field overcomes a critical threshold value. We will briefly discuss here the first 
approach, and refer to the next sections for the second one, which is the one 
taken in this is review. 

The attempt to model the Barkhausen phenomenon in terms of spin models 
is indeed appealing. Thanks to universality on expects the statistical properties 
of BN to be independent of the microscopic details of the system, provided that 
the relevant symmetries of the problem are properly taken into account. This 
allows in principle to derive macroscopic properties without going through the 
analysis of complicated micromagnetic equations, which are very difficult to 
treat analytically. Instead, one models the avalanche dynamics with simplified 
rules for the evolution of integer valued spins. This approach allows both to 
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perform large scale simulations, and to derive some theoretical result. 

In particular, among all the spin models, the non-equilibrium random field 
Ising model (RFIM), proposed by Sethna et al. [63] as a model for avalanche 
dynamics in BN, has been widely studied both analytically and numerically 
[1,64-69]. In this model integer valued spins Si = ±1 are assigned to each site 
of a d— dimensional lattice, and interact with the nearest neighbors through a 
ferromagnetic coupling J. The spins are coupled to the external field H, and 
to a quenched random field hi extracted from a Gaussian distribution with 
variance R, which is associated to each site of the lattice, and is supposed to 
mimic the presence of disorder. The Hamiltonian reads 



where the first sum is restricted to the nearest neighbors. In the non- 
equilibrium zero temperature dynamics proposed by Sethna and co-workers 
[63], at each time step the spins align with the local field: 



where, again, the first sum is restricted to the nearest neighbors of the spin 
i. Starting from a configuration with all the spins down, the field H is slowly 
ramped from — oo to +oo. A single spin flip can cause nearest neighboring 
spins to flip, and trigger an avalanche. In the quasistatic limit the field H{t) 
is kept constant during the avalanches. 

For small values of R, the random field is not strong enough to stop the 
avalanches, and the hysteresis loop has a squared shape, with a large jump at 
a critical field H c at which the big avalanche starts. When R is large, only small 
avalanches are observed: essentially every spin flips independently when H(t) 
is large enough to overcome the local pinning field hi, plus the comparably 
small energetic contribution from the coupling J with the neighboring spins. 
The two regimes are separated by a critical value R c of the variance of the 
random field, for which the avalanches are power-law distributed. 

This model predicts that for the critical value of the disorder R = R c , at 
the critical point, in d = 3, and in the quasi-static limit, the exponents in the 
power-law distributions of size and duration of the avalanches are given by 
r ~ 1.6, and a ~ 2.05. These values have been obtained both by large scale 
numerical simulations [67], and with an epsilon expansion around the mean 
field solution, that has been be worked out analytically [66]. 

The main problem with this approach is that RFIM predicts scaling only 
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Figure 17. Domain structure of _Fe2iCo64-Bi5 amorphous alloy observed by scanning electron 
microscopy. This is the typical structure observed in soft ferromagnetic materials. (Reprinted with 
permission from [56]). [56] Copyright 1998 by the American Physical Society. 

at a critical value of the disorder, while there is not reason to expect this 
condition to be satisfied in experiments. Another problem is that the results 
for the critical exponents in the non-equilibrium RFIM are only marginally 
compatible with the values observed in experiments, where the exponents are 
around r ~ 1.50, a ~ 2.0, or r ~ 1.27, a ~ 1.5 depending on the range of 
dipolar interactions. 

More importantly, even a rough investigation of the qualitative behavior of 
the non-equilibrium RFIM shows that the domain structure is very different 
from the one observed in real magnets (see figure [T7|) . While in the RFIM 
avalanches nucleate in different points of the system and grow with a fractal 
spatial structure, in real soft magnets domain walls are rather flat [2]. 

In the RFIM the magnetization proceeds simultaneously by two different 
mechanisms, nucleation and wall displacement, while in real soft magnets (and 
in the central part of the loop, where BN is measured) the wall movement is 
the only relevant mechanism of magnetization reversal: the dynamic evolves by 
the motion of existing boundaries, rather then by nucleating new domains. As 
a consequence, the RFIM, as it is, is inadequate as a model for BN in soft mag- 
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nets, while it could be an interesting starting point to describe magnetization 
dynamics in hard magnetic materials. 

The non-equilibrium RFIM neglects dipolar interactions, which originates, 
both from free dipoles which are present on the sample surface, and from 
those on the magnetic interface. Magnetostatic dipolar interactions due to the 
presence of "free magnetic charges" on the sample boundaries, generate a long 
range demagnetizing field, which opposes the external applied field. The most 
relevant ingredient missing in the RFIM is indeed the demagnetizing field, 
which reflects in the flat domain domain structure observed in soft materials 
(see figure [T7|) . very different from the structure found in the RFIM, where 
demagnetizing forces are not taken into account. The other crucial effect of 
the demagnetizing field, as we discuss below, is that it prevents the formation 
of large avalanches. 

Several variants of the RFIM have been proposed to overcome these prob- 
lems. A particularly successful one, the front propagation model (FPM), orig- 
inally introduced in [70], has been extensively analyzed in [51,69,71] as a 
model for BN. In this variant of the RFIM, the dynamics is such that only 
spins close to an existing domin boundary take part to the dynamics, while 
spins not adjacent to an avalanche front are prohibited to flip, even when it 
would be energetically favorable, so that the nucleation of new domains is to- 
tally suppressed. This can also be implemented by introducing in the RFIM 
Hamiltonian a long range term mimicking the demagnetizing field. The FPM 
gives exponents r ~ 1.27, and a ~ 1.72 in d = 3, which are compatible with 
BN experiments for soft magnets in the short range universality class. 

It has to be stressed that, unlike the RFIM, the FPM exhibits critical be- 
havior in the whole low disorder regime R < R c and not just at the critical 
value R c . Indeed, the demagnetizing field opposes the reversal process and 
forbids the magnetization to change abruptly, as it is the case, below R c , in 
the original non-equilibrium RFIM: the demagnetizing field reduces the effec- 
tive field until it stops the domain growth, which will start again with a new 
avalanche only by the increase of the external field. This mechanism will be 
treated extensively in the next section. The criticality observed is therefore 
not related to the disorder induced transition of the non-equilibrium RFIM, 
but to a depinning transition. 

Dipolar interactions between free dipoles on the magnetic wall produce an 
effective elasticity between different parts of the domain wall, that can be 
short range or long range, depending on the material. The long range elastic 
interaction has been claimed to lower the upper critical dimension to d=3 [56] . 
In such materials, a mean field approximation is therefore expected to apply. 
In mean field, the front propagation model gives r ~ 1.5, in good agreement 
with experiment on materials in the long range universality class. 

Another point that should be considered carefully in any RFIM approach 
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to BN regards the kind of disorder it represents, since random fields are not 
present in real magnets. It is therefore crucial to understand whether the RFIM 
non-equilibrium phase transition represents a broad enough universality class 
that is preserved when other kinds of disorder are introduced. This question is 
still controversial. Several studies have attempted to understand by numerical 
simulations whether other disordered spin models belong to the same univer- 
sality class as the RFIM. Examples include the non-equilibrium random bond 
Ising model (RBIM), with Hamiltonian 



where the disorder enters into the exchange couplings Jij, and the non- 
equilibrium site diluted Ising model (SDIM), described by 



where the quenched variables Cj, taking the values [0,1], indicate whether a 
given site represents a magnetic spin or not, thus modeling the presence of 
non-magnetic inclusions in the material. 

Simulations show that both the non-equilibrium RBIM [72] and the non- 
equilibrium SDIM [73] exhibit a disorder induced phase transition, with a 
behavior qualitatively similar to that arising in the RFIM. The critical expo- 
nents for the avalanche dynamics differ slightly, but are compatible, with those 
of the RFIM, suggesting that RFIM may indeed represent a broad universality 
class. The robustness of RFIM with respect to different kinds of disorder is also 
discussed by Dahmen and Sethna in [66], where it is argued that universality 
is expected to hold on the base of renormalization group calculations. The 
authors suggest that although RBIM and the random anisotropy model have 
higher symmetries then the RFIM in equilibrium, they have the same sym- 
metries in non-equilibrium conditions, since the external field at the critical 
point breaks the up-down symmetry. 

To conclude our short discussion of spin models, we can summarize that the 
simple RFIM is inadequate to describe the Barkhausen effect in soft magnets, 
rather it could be a good starting point to model magnetization reversal pro- 
cess in hard magnets. The FPM, obtained from the RFIM by incorporating 
the demagnetizing field, and suppressing domains nucleation, is instead a more 
adequate spin model for soft magnets, and give results in good agreement with 
experiments on materials in the short range universality class for d = 3, and 
for the long range in mean field. 

It is however very important to stress that in the FPM, the observed critical 
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behavior is not related to the proximity of the RFIM non-equilibrium critical 
point (indeed criticality is observed for all R < R c , i.e. in the whole low disorder 
regime), but rather to the depinning transition of the magnetic interface that 
will be discussed extensively in the next sections. 

In some sense spin models can be more general then interface models, as 
can easily deal for example with interface overhangs, however, in the case 
of soft magnets, given the nature of the transition, a simplified approach in 
terms of domain wall dynamics, seems more convenient, and has the advantage 
to allow to treat separately domain wall propagation and nucleation. The 
FPM naturally provides a discretized way to implement the front dynamics, 
very suitable for simulations. An alternative approach, which we follow in this 
review, is to formulate directly an equation of motion for the propagating wall. 



3 The ABBM model 

The first important theoretical approach to BN was done in 1990 by Alessan- 
dro, Beatrice, Bertotti, and Montorsi, who proposed a phenomenological 
model, which was then named after the authors of the two companion pa- 
pers [31,32], one dealing with theory, the other with experiments. By-passing 
all the difficulties that appear when trying to approach the problem from basic 
micromagnetic equations, Alessandro et al. managed to capture the essential 
ingredients at the origin of BN, by formulating a very simplified equation for 
the domain wall dynamics. Their model turned then out to be extremely suc- 
cessful in reproducing the experimental observations, and played a fundamen- 
tal role in most of the successive theoretical progress, and in the development 
of a consistent theory for BN. 

The ABBM model was inspired by the work of Neel [74, 75], who was the 
first to introduce a random energy model in the study of hysteresis: in Neel's 
model, the complex domain structure usually present in a magnetized sample 
is replaced by a single domain wall, moving in a random energy landscape, 
which is supposed to account for all the interactions among walls, and of the 
walls with the medium. The energy profile is constructed with a sequence 
of equispaced parabolic arcs, whose curvatures are extracted from a Gaussian 
distribution, in such a way that, on small scales, the motion is reversible, while, 
on large scales, the wall moves in random irreversible jumps. This model was 
introduced to explain in simple terms the Rayleigh law in hysteresis. 

In the Alessandro et al. approach, the idea of a random energy landscape 
is generalized, and used to construct a stochastic equation for the domain 
wall dynamics, in order to describe the Barkhausen effect. The equation is 
derived by following a previous work by Williams et al. [76], who studied the 
dynamics of a single domain wall under an external field in a FeSi single 
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magnetization 



Figure 18. Random energy landscape representing the pinning field, the arrows represent the 
Barkhausen jumps. In random energy models, when the system has enough energy to overcome a 

barrier, it jumps irreversibly to the next valley, with a corresponding random jump in the 
magnetization. Then the system is stuck until the effective field has grown enough to overcome the 

new energy barrier. 

crystal sample. In this work, the motion of the wall, limited by eddy currents 
damping, was described by a linear relation between the average velocity and 
the net field acting on the interface, given by the difference between the applied 
field H(t) and some threshold field Hq, which includes all kinds of random 
(like interaction with pinning centers) and non-random (like dipolar forces 
and demagnetizing fields) interactions (see figure [T8|) . 

The ABBM model is derived by making a more stringent claim, namely by 
assuming that the linear relation between velocity and effective fields holds 
not only on average, but also on the instantaneous velocity. By separating in 
Hq the random contribution of the pinning field H p from the demagnetizing 
component Hd, one could then write 



The further assumption is that the pinning field is Brownian correlated (see 
Eq. [T5l below), as suggested by some experimental evidence in special system 
containing just one active wall [77,78]. Let us point out that this assumption 
is quite crucial, since an uncorrelated pinning field, as the one considered 
in the original work by Neel, would not give rise to power law distributed 
magnetization events. 

In the following section we will first describe the ABBM equation, as it was 
formulated in [31,32]. Then, in order to put this model in a more general 
framework, we turn to the description of the dynamics of a magnetic elastic 
interface in a random media. Following [79], we describe a depinning interface 
model for BN, that in the mean field approximation gives back the ABBM 
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model [56]. As we will see, this will also clarify the origin of the correlated 
effective pinning field. 

3.1 ABBM as a phenomenological model 

The ABBM theory assumes that there is just one single rigid wall separating 
two domains of opposite magnetization. The system has therefore just one 
degree of freedom, and the dynamics is described in terms of a one-dimensional 
equation of motion for the center of mass of the wall. The wall moves under the 
action of an applied driving field in an effective Brownian correlated pinning 
field, mimicking the presence of disorder in the sample, as discussed above. 

Without loss of generality, we choose the origin of the coordinate indicating 
the position of the wall at the center of the sample, i.e. at the zero magne- 
tization domain configuration, so that, in what follows, the magnetization is 
proportional to the position of the wall. The equation of motion for the center 
of mass of the wall is then an equation for the magnetization. 

The motion is assumed to be overdamped, so that the velocity of the wall, 
which is proportional to the derivative of the magnetization, is equal to the 
sum of the forces acting on the interface, with a damping coefficient T, char- 
acterizing the amount of dissipation: 



As mentioned before, the total force acting on the wall, on the r.h.s. of the 
equation, includes two terms, besides the time-dependent applied magnetic 
field H(t), as in equation (fT3l) . One is the demagnetizing field Hd generated 
by free magnetic charges on the boundary of the sample, which opposes the 
external field and here is assumed to be proportional to the magnetization 
Hd = —km, through a constant k (demagnetizing factor) that takes into ac- 
count the geometry of the sample. The other one is the random pinning field 
Hp = W(m), which accounts for defects and non-magnetic impurities present 
in the sample, and it is assumed to be Brownian correlated in the magnetiza- 
tion: 



3.2 Limitations of the ABBM model 

The ABBM model has the advantage to be exactly solvable, and to repro- 
duce with striking accuracy most of the statistical properties observed in the 
experiments for soft magnetic materials in the long range universality class, 




(14) 
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described in section [2~4l However, it is worth to underline its main limitations. 
One, that we will deal with in the last part of the review (see section [5]), is that 
it is unable to reproduce some features of non-universal origin, in particular, 
the asymmetric shape of Barkhausen pulses. 

Another important limitation of the ABBM model is that it deals with a 
single domain wall, while in common experiments more interacting domains 
are present in the specimen and participate in the magnetization process. 
The ABBM predictions on the signal power spectra are for example cleanly 
observed in experiment where a single domain wall is present [78] , while other 
experiments suggest that the power spectra changes when more domains are 
present [80,81]. To fully understand these experiments a model which includes 
many interacting domain walls would be needed. This point was indeed raised 
in the original paper [31], where the authors argue that in the case of many 
interacting domains, the same equation of motion approximately holds, once 
the correlation associated to the pinning field (see the discussion in section 
I4.8P is multiplied by the number of active walls. 

Another feature of BN for which the description in terms of the ABBM 
model is unsatisfactory is the partial signal reproducibility observed in some 
experiments [82,83], that could be due to small variations in the initial condi- 
tions or in the driving field rate or in the peak field, or to thermal effects. A 
model able to describe this behavior should treat many interacting degrees of 
freedom explicitly, and not in an effective manner as in the ABBM model. 

A comparison point to point between ABBM model predictions and the ex- 
perimental evidences, which are summarized in subsection 12.4.91 is postponed 
to subsection 14.91 



3.3 Elastic interfaces in a random media 

An analysis of the Barkhausen phenomenon at a mesoscopic level is done by 
describing the system as an elastic interface in a random media driven by an 
external field, and relates the criticality observed in the experiments to the 
well known depinning transition of the interface. 

The response of an elastic manifold in a random media to an external driv- 
ing force is a well studied problem in statistical physics, that arises in many 
contexts, for example in the study of fluids in porous media [84-88], flux lines 
in type II superconductors [89], charge density waves [90-92], and crack fronts 
in solids [25,27,28,93-95]. In the following we briefly review the main facts, 
for the case of a d — 1 interface in d dimensions, and refer to the extended 
literature for more details [96,97]. 

Assuming that the motion is overdamped, so that inertial effects may be 
neglected, the general equation describing the dynamics of an elastic interface 
in a quenched random media, is obtained by equating the damping term to 
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the total force acting on the interface: 

where, assuming that there are no overhangs, z defines the position of the 
interface as a function of the d — 1 dimensional coordinate f on the substrate 
plane, and of time t. The coefficient T is the damping coefficient, ruling the 
rate of dissipation. 

The simplest equation of motion for an interface in an isotropic medium is 
the quenched Edward Wilkinson equation in a driving field, given by 

dz 

T-^(r,t) = F(t) + vV 2 z(r,t) + rj(z(r,t),r) , (17) 

where F is a uniform, eventually time dependent, applied force, V 2 z is the 
elastic term, whose coefficient v is the surface tension, and 7/(21, r) is the noise 
term, that mimics the presence of quenched disorder in the system. We con- 
sider here the zero temperature limit, where thermal noise is negligible with 
respect to quenched noise. A more general equation should include a kinemat- 
ically generated nonlinearity of the Kardar-Parisi-Zhang type, which changes 
the asymptotic scaling. However this term vanishes at the critical point and 
therefore Eq. [T7]is expected to have the correct critical behavior at the depin- 
ning transition (see below). 

The quenched disorder is described by a random potential V(z,f), whose 
derivatives gives the local pinning field r](z,f) in equation (|17p . In d = 3, if 
Pi = (xi,yi,Zi) are the coordinates of the point-like defects, and fj = (xi,yi), 
then the random force is 

r ] (z,r) = -RY,S(z-z i )5 2 (f-r i ), (18) 



where the sum runs over the pinning points. Here we assume the random noise 
to be uncorrelated and Gaussian distributed, so that 

(rj(z, r) V {z', ?)) = 25(z - z')5 2 (f- f*) (19) 

Systems described by equation (fTTl) exhibit a second order dynamic phase 
transition, called the depinning transition. The order parameter is the velocity 
of the interface, and the control parameter is the applied force, that we assume 
uniform and perpendicular to the substrate plane. 

The behavior of the system is sketched in figure [I9j At small applied forces 
the interface is locked by the pinning potential and v = 0. At some threshold 
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Figure 19. Sketch of the behavior of the interface velocity as an order parameter for the second 
order non-equilibrium depinning transition. 

value F c the force overcomes the pinning. For F > F c the interface moves with 
constant velocity. Just above the threshold value F c the interface moves in 
power law distributed jumps. 



3.4 Dynamics of a magnetic elastic interface 

In the case of a magnetic wall, the interface separates two domains with oppo- 
site magnetization. When an external field H is applied, the system tends to 
increase the domain with the largest component of the magnetization in the 
direction of the field, by domain wall displacement. Thus the field H acts as 
an effective force perpendicular to the wall, as shown schematically in figure 
1201 and we can replace F with H in equation (|17j) . 

Unlike many other cases where an elastic interface is driven in a disordered 
media, magnetic systems do not generally experience a depinning transition 
of the kind predicted by equation (fT71) and described in the previous section. 
This is due to the presence of other important interactions peculiar to the 
magnetic case, that modify the behavior of the interface in an essential way. 
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Figure 20. Magnetic wall in an external field H. The system tends to increase the size of the 
domain magnetized in the direction of the external field, which therefore acts as an effective force 

perpendicular to the interface. 



A first important step to adapt the general theory of elastic interfaces in 
random media in order to take into account the magnetic nature of the inter- 
actions was done by Urbach, Madison, and Markert [98] who modified equation 
(|17p to include the effect of the demagnetizing field (see also [99] ) . The demag- 
netizing field is the field generated by the presence of free magnetic charges on 
the boundaries of the sample (see figure [2lTa)). It opposes the applied field, 
and it is roughly proportional to the magnetization of the sample. The de- 
magnetizing field is assumed, as in the ABBM model, to be of the form —km, 
uniform and proportional to the magnetization, where the demagnetizing fac- 
tor k takes into account the geometry of the sample. From micro-magnetic 
calculations it can be shown that the assumption of a uniform demagnetizing 
field is strictly correct only in ellipsoidal geometries, however, it turns out to 
be a good approximation also in more general cases. 

The addition of the demagnetizing field to the equation of motion has a 
crucial effect on the dynamics, since it provides the restoring force needed 
to keep the interface constantly at criticality. Without this term the system 
would be critical just at a threshold value of the applied magnetic field, and 
experience a true depinning transition, while, in BN experiments, criticality is 
observed in an entire range of values of H(t) around the coercive field. These 
range of values correspond to the linear region of the hysteresis loop, which 
is precisely where the magnetization reversal is dominated by domain wall 
motion. 
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Figure 21. (a) Free magnetic charges on the boundary of the sample, (b) Free magnetic charges on 

the interface. 

The drastic effect of the demagnetizing field on the behavior of the system is 
confirmed by performing experiments in particular geometries, as Barkhausen 
measurements on toroidal samples. In this case closed interfaces can be forced 
so that, due to the absence of free magnetic charges on the boundaries, the 
demagnetizing factor is zero, and indeed, a true depinning transition is ob- 
served. This behavior has first been observed on a FeSi single crystal in a 
frame geometry [76], and later on several other materials [100-102]. In any 
other geometry, where k ^ 0, the equation of motion for the interface has to 
be modified with respect to equation (fT71) by replacing the applied force F(t) 
not just with the applied field H(t), but with an effective field 

H eff (t) = H(t) - km{t) . (20) 

The demagnetizing field produces the feedback mechanism that guarantees 
the presence of criticality without any tuning of the driving field. The mecha- 
nism is the following: the interface is stuck as long as H(t) is such that H e fj{t) 
is below the critical value H c . As soon as H e ft exceeds the critical value for the 
depinning transition, the interface starts to move in a series of discontinuous 
jumps. The displacement of the interface, however, causes an increase in the 
magnetization m(t), and therefore in the demagnetizing field, so that at some 
point the effective field is back below the threshold value, and the avalanche 
stops. 
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The work by Urbach et al. [98] represents a fundamental step towards the 
realistic description of an elastic magnetic interface, since the introduction of 
the demagnetizing field changes dramatically the behavior of the model, from 
one that shows criticality at the pinning threshold, to one that is critical in a 
finite range of values of the applied field. 

However, the equation used in [98] to describe the domain wall dynamic 
still does not take into account all magnetic interactions, namely, it neglects 
the possible presence of free magnetic dipoles on the wall (see figure [21] (b)). 
Indeed, the local curvature of the interface may give rise to discontinuities 
in the normal component of the magnetization, and generate dipolar forces 
originating in the bulk of the material. 

This complicated interaction has been worked out and included in the model 
by Cizeau, Zapperi, Durin, and Stanley [56,79]. Cizeau et al. considered the 
effect of dipolar forces originating from free magnetic charges on the wall, and 
evaluated this bulk effect, showing that it gives rise to a further contribution 
to the total applied force, acting as a long range elasticity, which works to 
flatten the interface. 

The long range elasticity term turns out to be of the form: 

J dnK{r- f{) (z(n) - z{f)) , (21) 

with a kernel K given by 

This interaction has been calculated under the mild assumption that the sur- 
face charges are not strong enough to deviate the magnetization from the easy 
axis. The coefficient Kq is given by Kq = ^Qm 2 /2iT, where n$ is the vacuum 
permeability, and m s is the saturation magnetization per unit volume. 

The full equation of motion, including the contributions from surface and 
bulk dipolar interactions finally results to be: 

Tz(r, t) = uV 2 z(f, t) + J cPnKtf-fi) (z(f{) - z(r))+H(t)-kz(t)+r](z,r) , 

(23) 

where z(t) denotes the average position of the interface and is proportional to 
the magnetization. 

The long range elasticity term is relevant in a renormalization group sense, 
since it is dominating the surface tension elastic term, as follows by power 
counting: rescaling x — > x' = bx and assuming that h — > h! = b a h with 



(22) 
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some positive a, the elastic term from surface tension is rescaled by a factor 
b a ~ 2 , while the long range elasticity term from dipolar interactions is rescaled 
by £> d+Q_3 , where d — 1 is the dimension of the substrate plane. This point 
is treated in a more rigorous way in [56] where the authors, following the 
functional renormalization group developed in [87,92,103], obtain the value 
d c = 2 + 1 for the upper critical dimension. The lowering of the upper critical 
dimension of the system from 4 + 1 (which is the one characterizing a generic 
elastic interface in a random media i.e. without the long range kernel) to 2 + 1 
as a consequence of the presence of the long range dipolar interactions ensures 
that the mean field approximation holds down to d = 3, and ultimately allows 
to describe most experimental situations in terms of a mean field theory. 

Corrections to scaling may however be present in d = 2 + 1: indeed the mean 
field power law behaviors are usually modified at the upper critical dimension 
by multiplicative logarithmic factors, as it is for example in the simplest case 
of the equilibrium Ising model in d = 4 [104]. As an example closer to the 
system described by Eq. [23"1 let us mention that logarithmic corrections at 
the upper critical dimension for the depinning transition have recently been 
calculated for the quenched Edward Wilkinson equation with a driving force 
(HID in [105]. 

Equation f|23|) contains all the energetic contributions involved in the domain 
wall dynamics. In principle, one should add an additional noise term to the 
equation, to account for thermal fluctuations. However, it has been observed 
experimentally [98] that temperature does not affect Barkhausen measure- 
ments in bulk materials. At this regard see also [56], where a rough estimate 
of the signal produced by thermal effects is given, which results to be well be- 
low the typical instrumental background noise. This could be different in thin 
films, where thermal activated motion has been observed. Recent experiments 
by Kim, Choe, and Shin [106, 107] report measures by magneto optical imag- 
ing on Co thin films, where the applied field is raised just below the coercive 
field and then kept constant, so that the observed Barkhausen jumps are only 
triggered by thermal noise (note that the Barkhausen jumps going outside 
the imaging window makes unreliable the estimates of the critical exponents 
obtained in [106,107], as already observed in [5]). therefore, in the case of thin 
films, an extra time dependent random term may be required in the equation 
of motion to include thermal noise. 

The elastic interface model described by equation ([23]) has been extensively 
studied by numerical simulations, by direct integration and by simulating an 
automaton version of the model, and both in finite dimensions and in the mean 
field limit. Simulations in the infinite range limit confirm the equivalence with 
the ABBM model (see below), by reproducing the same amplitude distribution, 
and pulses statistics. The simulations of the interface model in d = 3, recover 
the mean field results, in agreement with d c = 3 being the upper critical 
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dimension. However, the scaling with k of the cut-off in the distributions of 
pulses sizes and durations does not agree with the predictions of the mean 
field. This point will be discussed more in detail in subsection 14.101 

Although in this review we are dealing only with soft magnetic materials in 
the long range universality class, it is worth to note that the model by Cizeau 
et al. [79] has a more general validity. The results discussed above change when 
the long range kernel in equation f)23() is absent, and the surface term is the 
dominating term. In this case the model is equivalent to the one previously 
proposed in [98]. The upper critical dimension in this case is d c = 4 + 1, and 
simulations in d = 3 give r ~ 1.3 and a ~ 1.45 for the critical exponents char- 
acterizing the pulses' distributions [108-110]. Simulations results are confirmed 
by renormalization group calculations [87, 108, 111] which give r ~ 1.25 and 
a ~ 1.43. In this case, no dependence on the driving field rate c is observed. 
These results are in very good agreement with experiments on materials in the 
short range universality class [56,112], indicating that the magnetic elastic in- 
terface model may indeed have a very general validity, providing a description 
for the magnetization reversal process suitable to deal with both universality 
classes of soft ferromagnets. 

3.5 The ABBM model as a mean field theory 

A mean field version of the model by Cizeau et al. [79] has been worked out 
by the same authors in [56], where they show that in the mean field approx- 
imation, their interface model coincides with the phenomenological ABBM 
model. This fact might indeed explain why the ABBM model describes the 
experiments so accurately. 

Mean field theory is derived by taking the limit of long range interactions in 
the equation of motion. This is done by coupling all the points on the interface 
uniformly with the average wall position. Equation (|23|) become 



where the coefficient J includes the contributions from both short and long 
term elasticity. By averaging this equation over r, one gets an equation for the 
average position of the interface: 



where S is the surface of the substrate plane. The integral over the noise gives 
rise to an effective pinning W(z), which depends only on the average position 
of the interface. 




(24) 




(25) 
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Indeed from equation (fT8j) 

,( 2 (f, ( ), r 1 = I i^mU^ = -fl^- W , t )-„), (26) 

J i 

which is proportional to the number of defects points that are pinning the 
interface to the average position ~z. The effective noise W is then given by 

W(z) = ± j d 2 rr ? (z(r,i),f) = -|^,5(z(r,t)-z,)- (27) 

i 

To calculate the correlations, note that, in the difference W(z) — W(z'), the 
only terms that contribute to the sum are those relative to the pinning points 
that have actually moved, whose number is proportional to the average wall 
displacement z — ~z', so that, r\ being random and uncorrelated, the effective 
noise is a Brownian noise in z: 

((W(z)-W{z J )) 2 } = D{z-z'). (28) 

This finally explains the nature of the effective Brownian correlated pinning 
field, which was observed experimentally in [77,78], and used as an assumption 
in [31] to derive the ABBM equation. The long range correlations observed in 
the pinning field are not a result of real long range correlations between defects 
in the material. Rather, the effective correlated pinning field experienced by the 
center of mass of the wall originates as a collective effect from the interactions 
along the elastic wall moving in an uncorrelated disordered medium. 

The average wall position is proportional to the magnetization, so that equa- 
tion [25] can be rewritten as 

IYh(i) = H(t) - km(t) + W{m) , (29) 

which coincides with the ABBM equation (|14p . 

The ABBM model is extremely successful in describing BN experiments, 
although it was proposed on phenomenological basis, without a clear reference 
to the underlying physics. The fact that the ABBM model can be derived as 
a mean field theory of an interface model obtained by taking into account all 
the relevant magnetic interactions is therefore very important. We want to 
stress that the relation with the elastic interface model not only validates the 
ABBM model, by giving to it a more grounded microscopical derivation, but 
also allows to relate the model parameters to microstructural properties of the 
material, making possible a close quantitative comparison between theory and 
experiments. 
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4 Mapping ABBM onto a random walk process: exact results 

The ABBM approach allows to write a simple equation of motion for the center 
of mass of the magnetic wall separating two regions of opposite magnetization. 
Interestingly, the statistical properties of this dynamics can be fully worked 
out analytically. 

The easiest way to do this is to transform the equation by eliminating the 
correlated noise. We show below that equation (|14p maps onto an standard 
Langevin equation for a Brownian motion in a suitable potential. The equation 
has to be solved taking care of the boundary conditions, that ensure that the 
signal always stays positive. 

4.1 From the ABBM equation to the Langevin equation for a biased 
random walk in a logarithmic potential 

Equation (|29p takes a particularly simple form in terms of velocity versus 
magnetization. Taking a derivative with respect to time, and defining v = 
dm/dt 



where we used d/dt = vd/dm. 

Assuming that the field H(t) is varied at constant rate c, as in most BN 
experiments, and defining w(m) = dW/dm, we get 



where the noise term, being the derivative of a Brownian noise, is delta corre- 
lated in magnetization, with {w(m)w(m')) = 2D5(m — m'). 

Note that, the ABBM is defined with the assumption that just one single 
and rigid wall is present. This allows the identification of the magnetization 
with the coordinate defining the wall position, and therefore the variable v, 
which is the derivative of the magnetization with respect to time, can be 
interpreted both as the average velocity of the interface, and as the induced 
voltage directly measured in inductive BN experiments. 

The four parameters in equation ()29[) can be reduced to two by rescaling 
the velocity v = vT/D, and defining k = k/D, c = cT/D 2 , and r\ = w/D, so 
that {n(m)r](m')) = 25(m — m'): 



dv _ H(t) 




dW 



(30) 



dm v 




(31) 



dv c ~ , . 

- — = k + mm) 

dm v 



(32) 
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This leaves two independent parameters k, and c. In order to avoid cumber- 
some notations we will avoid the tildes in the following. 

Interpreting the velocity v as a coordinate x and the magnetization as time 
t, this is the Langevin equation for a random walk in a logarithmic plus linear 
potential U(x): 

U(x) = — clog x + kx , (33) 

where the coefficients of the two terms in the potential, that rules the dy- 
namics, correspond to rescaled driving field rate and demagnetizing coefficient 
respectively. 

The Fokker-Planck equation corresponding to the Langevin equation 

x = - - k + r)(t) . (34) 

x 



IS 



dt dx \ dx V x 

where P(x, t) is the probability for the walker to be at x at time t. 

The first thing that can easily be done is to solve equation [35] in the station- 
ary t — > oo limit. This gives 

Poo(x) oc x c exp(— kx) . (36) 

Expression [36] gives the time independent probability distribution for the 
process for all values of c, however, it will become clear further on that, once 
the boundary conditions are taken into account, the signal distribution can be 
obtained by simply taking the asymptotic solution, as has been done here, only 
in the high driving rate regime (large c) . We will come back to this distribution 
in the next section. 



4.2 The Fokker-Planck equation for a random walk in a logarithmic 
potential confined to the semiplane 

The linear term in the potential (|33p acts as a constant bias to the random 
walk, that, as we will see, has the effect to introduce a cut-off in the return 
times distribution. To proceed, we forget for now about the bias term, and 
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focus on the behavior of an unbiased random walk in a logarithmic potential: 

x = - + r)(t) . (37) 

X 

Interestingly, the same equation (|37|) in two dimensions appears in the anal- 
ysis of the annihilation dynamics of vortex-antivortex pair in the two dimen- 
sional XY model, and has been studied in that context by Bray in [33], which 
we follow in the next sections for what concerns the calculation of the persis- 
tence properties of the process. 

The Fokker-Planck equation associated to the Langevin equation (|37ft is: 



dP(x,t) = d_ f dP(x,t) _ c \ 
dt dx\dx x ' / 



(38) 



where P(x,t) is the probability for the walker to be at x at time t. 

Since by definition the avalanche ends when the velocity is back to zero, 
we are interested to the solutions of this equation with absorbing boundary 
conditions at the origin x = 0: 

P(0,t) = 0. (39) 

Given the absorbing condition, a walk starting at x = would never leave 
the origin, thus we choose as initial condition for the walk a small e << 1: 

P(x,0) = 5(x-e) (40) 

and take the limit of vanishing e at the end. 

The solution of equation ([38]) can be searched by separation of variables [33] 
in the form 

P(x, t) = r(x) exp(-k 2 t) , (41) 

where r(x) has to satisfy 

x 2 r"(x) — cxr'{x) + (c — k 2 x 2 )r(x) = . (42) 

Equation (|42p can be reduced to the Bessel equation by the change of variable 
r(x) = x a R{x): 



x 2 R"(x) + (2a - c)xR'(x) + ((a - l)(o - c) + k 2 x 2 )R{x) = , (43) 
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by choosing for a the value a = (1 + c)/2: 

x 2 R"{x) + xR'(x) + (k 2 x 2 - v 2 )R{x) = , (44) 

with v = ((1- c)/2). 

The general solution of equation (|38p can therefore be written as a linear 
combination of Bessel functions as: 



f' c 

P(x,t)=x l - U 

Jo 



dk (a(k)J u [kx] + 0(k)J- v [kx]) exp(-k 2 t) . (45) 



The coefficients a(k) and (3{k) are fixed by imposing the initial condition 
equation (|40p . that can be enforced by using the orthogonality properties of 
the Bessel functions: 



CO 



dkkJ a [kxi] J a [kx2\ = S(x± — X2)/x\ . (46) 



o 



This gives 



/ S \ i-f r°° 
P(x,t | e,0) = e f-J / cffcfc^J^eJJ^fcr] + £ J-^fce] Jl^kc]) exp(-/c 2 t) 

where A, and B are numerical coefficients. Finally, by performing the integrals 



e fx 



l-v ( x 2 + e 2 



t\eJ \ it 



P(x,t\e,0) = ~[~) exp \[AI V - +BI- V - ) , (48) 



ex 



2t 



ex 



2t 



where I u [z] is the modified Bessel function. 

To simplify the notations let us rewrite equation (|48p as 



P(x, t | e, 0) = AP^z, t | e, 0) + BP- v (x, t\e,0) (49) 

and analyze the two terms separately. 
For what follows is useful to introduce the probability current 

j(x,t) = -(d x P-cP/x) (50) 

so that dtP = —d x j. The absorbing boundary condition equation ()39l) corre- 
sponds to requiring a negative current at the origin. 

In order to evaluate the probability current at the origin corresponding to 
the two terms in equation (149 j) . it is convenient to expand the Bessel functions 
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in the limit of small e, using 

CO 

I u [z] = {z/2f J2 (z/2) 2m /(mT[m + v + 1]) . (51) 

m=0 

The first term in equation (|48p simplifies to 

P»(x,t | 6,0) = f ^-l^(4t)-^) e 2 ^exp (-^j (l + 0(e 2 )). (52) 
The corresponding probability current is 

3u(x,t | 6,0) = f ^i_ I (4t)-( 1 +") e ^ ((c- 1) + ^) exp (-^j (l + 0(e 2 )) , 

x / (53) 

and gives a non-vanishing contribution at the origin, which is negative, as 
requested for an absorbing boundary, as long as c < 1. 
The second term in equation (|48p simplifies to 

P-*(x, t | 6, 0) = ^L-^AtT^x 1 - 2 " exp (-^j (1 + 0(e 2 )) . (54) 

The coefficient 1 — 2v — c of the leading term of order x~ 2u in the probability 
current vanishes, and the next-to-leading term gives: 

4 ,r 2 " 2iy / r 2 \ 

j_ u (x,t | 6,0) = -^-^-a-") — exp (- Tt ) (1 + 0(6 2 )). (55) 

Since c is a positive quantity, the exponent 2 — 2u = 1 + c is always positive in 
the region of parameters of physical interest, and the contribution of the term 
P- v to the probability current at the origin vanishes. Therefore the physically 
relevant solution to our problem, satisfying the absorbing boundary condition 
is P u , and we take P = P v and drop the u sign from now on. 
Finally, we can conclude that the total probability current at x = is: 

j(0 , t ) = - l^l^L (4t)-( 3 - c )/ 2 6 1 - c , (56) 
j\ , ) r[(3-c)/2] v 1 ' y ' 

where we have replaced v with its expression as a function of c. The current 
goes to zero as c approaches 1, correctly suggesting that, as we will see more 
clearly in section [475], c = 1 separates two different regimes. 
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It must be stressed that for c < 1 the probability P(x, t | e, 0) in equation 
([52]) is not normalized to unity, when integrated over x. Indeed, due to the 
absorbing boundary, part of the probability density flows to the origin, so that 
the conserved quantity is 



where P a bs{t I e > 0) is the probability that the random walk has been absorbed 
at some time t' < t. This probability goes to 1 as t — > oo. For this reason, the 
overall time independent process amplitude distribution cannot be obtained 
by taking the t — > oo in P(x,t | e, 0), since the asymptotic distribution is zero 
for any x ^ 0. 

The asymptotic distribution clearly is not the correct quantity to be com- 
pared with the signal amplitude distribution in equation ([5|). A more mean- 
ingful probability, which indeed corresponds to the amplitude measured in 
experiments is obtained by integrating P(x, t | e, 0) over t. In general, for 
a stationary process, the asymptotic distribution P 00 (x) coincides with the 
average over time of the probability P(x,t): 



However, the quantity in equation ([5]) that is measured in experiments is the 
average over many avalanches of the probability for the signal to have a given 
value. During the avalanche, the process described by equation [3H is non- 
stationary, and thus the integral over t cannot be replaced by the asymptotic 
distribution. 

This point needs a further remark. In subsection 12.31 we stressed the fact 
that the BN experiments has to be restricted to the region where the signal is 
stationary, and this may generate confusion when compared with the present 
argument. The BN process is described by the stochastic equation ([30]) only 
when v > 0, that is as long as the effective field (applied field diminished 
by the demagnetizing field) overcomes the pinning field. As long as v > 
the magnetization increases, causing a variation in the effective field, which 
decreases until it goes somewhere below the critical depinning threshold. At 
this point, the interface stays locked, until the applied field brings the effec- 
tive field again above threshold. Therefore while v = 0, the process is just 
deterministic. Given that the wall is at rest, there is no noise term involved. 
The stationarity of the whole BN process (the random process at v > plus 
the deterministic one for v = 0) should not be confused with the stationarity 
of the process described by equation [34j which, in the low driving field rate 




(57) 




(58) 
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regime, is reached only when the random walk is absorbed by the boundary. 

Let us now go back to our problem of evaluating the signal amplitude dis- 
tribution. By taking integral over t in equation (|52[) to the first order in e we 
can extract the small x power law behavior: 



This expression is not valid for large x, where the small e expansion does not 
hold. The power law dependence on x of the process amplitude coincides with 
that of equation (f36|) obtained in the previous section, apart from the expo- 
nential cut-off depending on the bias term, and which has been neglected here. 
Indeed equation (|36|) gives an expression for the process amplitude distribution 
correct for all values of c. 

The large t asymptotic has been used in literature [109] to explain the ex- 
perimental power law in equation ([5]), however, we stress that this behavior 
cannot be obtained in general by simply taking the asymptotic solution. 

To compare this result with the experimental one on the distribution of the 
Barkhausen signal we must remember that here t corresponds to the physical 
variable m, and x to the wall velocity, or BN amplitude v. Going back to 
physical time will therefore produce an extra v _1 factor, recovering the c — 1 
exponent found experimentally (see equation ([5])). 

4.3 Distributions of sizes and durations of pulses from persistence 
properties 

Let us now turn to the analysis of the statistical properties of the signal pulses. 
The distribution of the sizes and durations of Barkhausen pulses can be derived 
from the persistence properties of the process (j34"|) . 

Let us introduce the first passage time distribution Po(0,t | e, 0), i.e. the 
probability that a random walk starting at some threshold value e at time 
t = goes back to e for the first time at time t, and the persistence probability 
P±(t | e,0), defined as the probability that a random walk that starts at e at 
time t = stays positive up to time t. 

These two quantities are related by the following relation: 




.2u+l 



= X 



,c 



(59) 




(60) 



since the probability for a random walk to be positive up to time t is equal to 
the probability that it comes back to for the first time at some later time 
t' > t, which is given by the r.h.s. of equation (|60p . 
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The persistence probability P±(t | e, 0) is also related to the solution of the 
Fokker Planck equation with absorbing boundary condition by 

POO 

P 1 (t\e,0)= dxP(x,t | e, 0) , (61) 
J o 

since, to survive up to some time t, a walker must be found at time t at some 
positive value x. This is equivalent to the normalization condition (|57p where 
P abs (t\e,0) = l-P 1 (t\e,0). 
Using the previous two equations we can calculate Po as follows: 

POO 

P (i,| e,0) = -dtP^t | e,0) = - / dxd t P{x,t | e, 0) , (62) 

Jo 

where we passed the derivative under the integral. 

Rather than replacing P(x,t \ e, 0) with its explicit expression and then 
dealing with the complicated integral, it is convenient to replace the integrand 
on the r.h.s. with its expression given by the Fokker-Planck equation: 

P (t,\ e,0) = -J dxd x (d x P(x, t\e,0)- ~P(x, t | e, 0)) = j(0, t) (63) 

where, in replacing the integral with the current in zero we use the fact that 
P(x,t | e, 0) and its derivatives decay exponentially to zero for x — > oo. Eq. 
(I63|) is consistent with the probability P aos for the process to having being 
absorbed before time t being equal to the total current flow through the origin 
up to time t: P abs = - J j(0, s)ds. 
Finally we get from equation (|56|) : 

P (t,\ e,0) oc (1 - c)t-( 3 ~ c )/ 2 . (64) 

When c — > equation [M] gives back the well known result for the persistence 
probability of a free random walk. 

To see how this relates to the distribution of sizes and durations of BN pulses 
one has to go back to the physical variables, where t has to be interpreted as 
magnetization. Equation (|64p gives the probability that the avalanche ends at 
some given value of the magnetization. It gives therefore directly the distri- 
bution of sizes of the avalanches, which is p(S) ~ 5 _ ( 3_c )/ 2 . This gives the 
correct exponent r = 3/2 in the quasi-static limit c = 0, and the expected 
linear dependence on c for finite driving field rates. 

The distribution of durations of BN pulses can be also derived by this same 
result, by imposing that p(S)dS = p(T)dT, and using the scaling S ~ T 2 
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that we will show in subsection 14.71 This gives p(t) ~ T~( 2 ~ c \ with a = 2 
for vanishing field rate, and a linear dependence on c that agrees with the 
fit of experimental histograms. Note that the exponential cut-offs cannot be 
recovered from this analysis, since they involve the demagnetizing factor that 
has been put to zero here (see beginning of section 14. 2p . 

4.4 Non-universality of the exponents and marginality of the 
logarithmic potential 

The continuous dependence on a parameter in otherwise universal exponents 
is usually related to the presence of some potential which is marginal in the 
renormalization group sense. Indeed, also in the present case, the dependence 
of the exponents a and r on the driving field rate c can be attributed to some 
marginal contribution. 

This becomes clear in the mapping onto the random walk process. The 
non-universality of the exponents in the persistence and first passage time 
probabilities are special to the logarithmic potential, and are related to the 
marginality of the 1 jx perturbation to a free random walk in a renormalization 
group sense. Indeed, analyzing the relevance of a power law perturbation x a 
to a free random walk, one can show by simple scaling arguments that the 
force is irrelevant for a > — 1 (the random walk behaves essentially as if it 
where free), while it is relevant, and dominates the noise term for a < — 1 
(the motion is in many respects deterministic, and persistence properties are 
governed by rare fluctuations of the noise) . 

The equation for a random walk in a generic power law potential is 

x = cx~^ + r](t) . (65) 

Under rescaling x — > ax, t — ► a z t, the coefficients rescale as c — > a z ~ 1 ~^c, 
D — > a z ~ 2 D, where D is the noise correlation. The equation of motion therefore 
becomes 

x = Co*- 1 -**-* + a z / 2 -yt) . (66) 

This equation has two fixed points, depending on the value of ip. The poten- 
tial term dominates with respect to the noise term for z — 1 — ip > z/2 — 1. In 
this case, the noise term renormalizes to zero, and the value of the dynamical 
exponent is obtained by balancing the potential term with the l.h.s., which 
gives z = 1 + ip. Replacing the fixed point value of z in the inequality we 
deduce that this is the fixed points for ip < 1. 

The other regime is the one where the noise dominates over the potential, for 
z/2 — 1 > z — 1 — tp. In this case, the coefficient of the potential renormalizes 
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to zero, and the value of the dynamical exponent is obtained by balancing 
the noise with the l.h.s., which gives z = 2. Replacing the value of z in the 
inequality we verify that the system renormalize to this fixed point as long as 
a > 1. 

Therefore, for tp > 1, z = 2, and c rescales to zero: the force is irrelevant, and 
the random walk behave essentially as if it were free. For ip < 1, z = l+^>, the 
noise correlation D rescales to zero: the noise term is (dangerously) irrelevant, 
and the process is almost deterministic. 

The case of our equation, tp = 1, corresponding to the logarithmic potential, 
is marginal: it is the value of ip for which all the three terms in the equation 
balance, with z = 2. This is at the origin of the non-universality of the per- 
sistence exponents and of their continuous dependence on the coefficient c of 
the potential. 

The effect of the driving field rate on crackling events has been recently stud- 
ied by White and Dahmen in [55] . In this work the authors make a distinction 
between an avalanche (sequence of events triggered by an initial event at con- 
stant driving field), and a pulse (uninterrupted sequence of events measured 
in experiments). In general, at non-zero sweep rate a pulse can be a superpo- 
sition of several simultaneously propagating avalanches, while avalanches and 
pulses coincide only in the adiabatic limit. Under some general assumptions 
that are fulfilled by many crackling systems, the authors of [55] show by scal- 
ing arguments that the value of the exponent a in the distribution of events 
durations determines whether or not the pulse statistic will be affected by the 
sweep rate c. In particular, they claim that no dependence on c has to be ex- 
pected for a < 2, while a linear dependence of the kind a = 2 — c is predicted 
in the case a = 2. It would be very interesting to make a connection between 
the results obtained in [55] in terms of superposition of adiabatic avalanches 
with those reported here for the ABBM model, but this has not been done 
yet. 



4.5 Existence of a threshold in the driving field rate from recurrence 
properties of the free random walk 

Another evidence in Barkhausen noise phenomenology is the fact that the in- 
termittent behavior typical of the avalanche dynamics is lost when the driving 
field is cycled too fast: above a threshold value of the driving field rate the typ- 
ical intermittent signal disappears, and the magnetization reversal proceeds in 
a unique avalanche. 

The existence of a threshold value in the driving field rate c, corresponds to 
the existence of a special value c for the coefficient of the logarithmic potential, 
that discriminates among two qualitatively different behaviors. The value c = 1 
already appeared in equation (I56j) as the threshold below which the probability 
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current at the origin is negative. 

Essentially, this value of c determines whether or not the random walk is 
recurrent. The random walk is certainly recurrent for c < 0, which corresponds 
to an attractive potential at the origin. This, however, does not correspond to 
the physical case of BN. The case of physical interest is c > 0, corresponding 
to a repulsive potential. Interestingly, the random walk may come back to 
the origin even when it is subject to a logarithmic repulsion, as long as the 
repulsive potential coefficient is not too strong. This gives a condition on the 
coefficient c. 

A a simple and clear way to interpret this behavior is to observe that there 
exists a natural mapping between the free random walk in d dimension and 
the confined one dimensional random walk in a logarithmic potential. Starting 
from the Fokker-Planck equation for a free random walk in d dimensions. 



we can write the Fokker-Planck equation for the radial part in polar coordi- 
nates as follows 



This equation shows that the radial part of free random walk in d dimensions 
performs a one dimensional random walk in a logarithmic potential confined 
to a semiplane, with coefficient c equal to d — 1 for the potential. 

Using this mapping, the existence of a threshold value in the coefficient of 
the potential is explained in terms of the well known recurrence properties of 
the free random walk. The condition d < 2 for the random walk to be recurrent 
translates into c < 1. 

4.6 Cut off in the distributions of pulses' sizes and durations as an 
effect of the bias to the random walk 

Let us now analyze the effect of to the bias term k. We consider a simple 
biased random walk first: 



By rescaling x — ► x + kt, it is clear that the probability for a biased random 
walk with constant bias k to return to the origin is equivalent to the probability 
for a free random walk to return onto the line kt. Given that the random walk is 
sub-linear, this introduces a characteristic timescale r, after which the random 



dtP(x,t) = V x P(x,t) 



(67) 




(68) 



x(t) = -k + r)(t). 



(69) 
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(a) 



(b) 
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Figure 22. Sketch of the mapping of the return of the biased random walk to the origin (a) onto 
the return of the free random walk to the line kt (b). The characteristic time r is denned as the 
typical time at which the free random walk crosses the line kt. 

walk will very unlikely, namely with exponentially small probability, go back 
to the line kt (see figure [22]) . 

The value of the characteristic time r is such that kr ~ y/r, which gives 
r ~ k~ 2 . Since t corresponds to magnetization in the physical variables, r 
corresponds to the characteristic avalanche size So in the cut-off of equation 

([6]). As we will see further on, To ~ S^ 2 , gi ym g ^o ~ k~ l for the cut-off in 
equation ((7|). 

As we have seen in subsection 14.41 the logarithmic potential term does not 
alter the dynamic exponent z which is still equal to 2 as for the free random 
walk, thus the walk in the logarithmic potential is still sublinear, and is affected 
by the bias in the same way as the free random walk. 

The fact that the bias term k generates an exponential cut-off in the dis- 
tribution of the return times, physically means that the size of the largest 
magnetization event is regulated by the demagnetizing field. 

4.7 Average shape of Barkhausen pulses as excursion of the random 
walk process 

In analogy to conventional critical phenomena one might expect that pulses 
of different durations have an average shape that, once properly rescaled, col- 
lapses onto a universal scaling function. In a recent paper by Sethna et al. [1] it 
was suggested that this scaling function is an interesting quantity to measure 
in order to make quantitative comparison between the theoretical models and 
experiments. This quantity provides a much sharper tool to test theory against 
experiments than the simple comparison of critical exponents, that are scalar 
quantities. 
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Indeed, by performing this kind of analysis on experimental data, average 
pulses, properly rescaled by their durations, approximately collapse onto the 
same function (see figure [Till [113]. 

The average shape of the pulse is related, in terms of the corresponding 
stochastic process, to what is known in the mathematical literature as the 
excursion, namely the average trajectory between two successive return of the 
process to some reference value. 

The average shape for the process (|37p has been calculated in [113]. This 
work was exactly motivated by the problem, raised in [1], of finding the average 
form of a Barkhausen pulse as a function of time. 

Let us consider a generic random walk in a potential U{x): 

±{t) = U{x)+r,{t). (70) 

The average over trajectories starting at x = at time t = 0, and constrained 
to come back to x = for the first time at t = T, that we will denote as {x(t))T 
is found to scale, in appropriate time regimes depending on the potential U, 
as 

(x(t)> T = ry{t/T) (7i) 

where / is some scaling function [114,115]. 

Let us consider processes constrained to the semiplane x > 0, and focus on 
the returns to the origin (take as reference value a small threshold x = e). 
Once the solution with absorbing boundary conditions of the Fokker Planck 
equation corresponding to the Langevin equation (J7DJ) is known, the excursion 
can be expressed by averaging over all trajectories going through x at time t. 

In general, the probability for a Markovian process starting at e at time 
and going back to e at time T to be found at position x at time t, is given 
by c(x,t;e, T) = P(x, t; e, 0)P(e, T; x, t). Using time translational invariance 
c(x, t; e, T) = P(x, t; e, 0)P(e, T — t;x, 0) = P(x, t; e, 0)P(x, T — t;e, 0), where 
P is the conditional probability for the inverse process. 

The average excursion is therefore obtained by taking the following limit for 
e -> 0: 

= ^ d xxP { x, t] e,0)P { x,T- t -,e,0) 

J™dxP{x,t;e,0)P(x,T -t;e,0) V ; 

which counts the fraction of excursions that go through x at time t. 

It is instructive to calculate the excursion for the simplest case of a free 
random walk [116]. The probability P(x, t; e,0) for the walk confined to the 
semiplane can be calculated by the image method as a linear combination of 
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two solutions of the Fokker Planck for the free process: 

P ( X) t ; e,Q) = J^ (exp-(— ) 2 /2* - exp"^) 2 / 2 ') , (73) 
which for small e gives 

P(x, t; e, 0) = J^-r 3 / 2 ex exp" x ' 2 / 2t . (74) 

In this case P = P due to the process time reversal invariance, therefore the 
average excursion is given by 

_ iT^exp-* 2 / 2 " 

\ X )T — r oo , o _ r 2 io T i K {t) ) 

J axx z exp x / ZT 
with t = t(T — t)/t. Evaluating the integrals gives 



Therefore the exponent 7 is equal to the wandering exponent of the free 
process 7 = 1/2, and the scaling function is 



f{x) = y/8fry/x(l-x), (77) 

normalized to have L dxf(x) = 1. 

Let us now consider the case interesting to our theory of the excursion for 
the random walk in a logarithmic potential. This can be calculated in a similar 
way, by observing that the process is still time translational invariant, and it 
is invariant under t — > —t, as long as the coefficient c is also inverted in the 
revers process. Therefore P(x,T — t;e, 0) is obtained by replacing c with — c 
in P(x,T-t;e,Q). 

Indeed, by repeating the previous steps with the solution of the Fokker Plank 
equation (f38l) . the calculation finally reduces, also in this case, to evaluating 
the ratio between two successive moments of a half Gaussian with variance 
w = 2t(T — t)/T, so that the average excursion results to be proportional to 
^/w = \/t(T — t)/T, as in the case of the free random walk, and it is still given 
by equation ([77]) . 

To translate this result into the shape of the Barkhausen pulse, we have 
to go back to our original physical variables. The excursion calculated above 
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Figure 23. Schematic representation of the average excursion. Different realizations of the 
stochastic process all returning at the origin after the same time interval T\ are averaged to give a 
smooth curve. Repeating the same process with a different return time interval T2 gives another 
curve, that can eventually be rescaled onto the first. 

corresponds to the signal v as a function of the avalanche size s. The average 
pulse as a function of s and of the total size S of the pulse is then 

(v(s,S))=S 1 / 2 f(s/S). (78) 

However, the commonly measured experimental pulses, as those shown in fig- 
ure Q31 represent the signal v as a function of duration t of the avalanche, that 
will have some different scaling and shape 

(v(t,T))=TPg(t/T). (79) 

Following [113] we calculate the scaling function g: we have to express s and 
S as functions of t, and T. By definition, the size s of the avalanche at time t, 
is the integral of the signal v from time to t: 

(s(t, T)) = ( f dt'v(t>, D) = T^ +1 dxg{x) , (80) 
Jo Jo 
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therefore, the total size S scales as S — T@ +1 . Approximating s(t,T) with its 
average value in the equation 

v(t,T)=v(s = s(t,T),S = S(T)) (81) 

gives the following condition on g: 

T?g{t/T) = T^' 2 f (J** dxgix)^ . (82) 

This enforces /3 = 1, and gives an implicit integral equation for g: 

g{x) = f ( f X g(x')dx') , (83) 



o 



that has to be solved subject to the boundary conditions g(0) = g(l) = 0. 
Equation (|83H takes a simpler form in terms of h(x) = Jq g(x')dx': 



h'(x) = f (h(x)) = iry/h(x)(l - h(x)) , (84) 

which is solved by h(x) = sin 2 (-7rx/2). This gives for g(x) 

g(x) = h'(x) = Tr/2sin(vrx) , (85) 

which is normalized to unit and satisfies ^(0) = g(l) = 0. Finally, for the 
normalized average avalanche as a function of its duration we get: 

( W (t,T)> = ^ sin rrf/T. (86) 

Note that, while Eq. [78] with / given by Eq. [77] is an exact result, and gives 
the scaling form of an avalanche as a function of magnetization, to derive Eq. 
I8B1 one has to replace the dependence of s on time in a single realization with 
its average, given by Eq. [80] 

Figure [Ml shows the comparison between rescaled experimental pulses of 
figure E] and the scaling function (|85jl predicted by the ABBM. The curves 
have similar shapes, however, the experimental collapse is characterized by 
an evident leftward asymmetry, while the ABBM predicts a symmetric shape. 
Note that normalizing the average pulses with iV as done in figure l24l should 
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Figure 24. The approximate collapse of average pulses in figure [14] is compared with the prediction 
of the ABBM model (bold line) from equation 1851 (Adapted from [5]). (Reprinted with permission 
from [5]). [5] Copyright 2005 by Academic Press. 

in principle be equivalent to rescaling them by T _1 , as in equation (|86p . since 



and S ~ T 2 . However, if one directly plots (v(t,T))T _1 , the curve collapse is 
obtained only for short avalanches, while a clear deviation appears for longer 
avalanches (see figure [25]) 

This implies that the scaling S ~ T 2 does not hold for large avalanches. 
Indeed, by plotting the average size (S) versus T for the same BN signal, the 
curve deviates from the predicted scaling for avalanches longer then ~ 0.6ms, 
which is more or less where the pulses start to deviate form the expected 
collapse (see figure [26]) . 

4.8 Power spectra 

The power spectrum for the ABBM equation has been calculated in the orig- 
inal paper [31]. In order to get the correlation function it is more convenient 
to use directly the physical variables, time, magnetization, and velocity (or 
signal amplitude). The probability R to find a signal v at time t is related to 
the probability P used in the previous sections by 




(87) 
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Figure 25. Average avalanche shapes for a polycrystalline FeSi material. Signals as a function of 
time v(t,T) are scaled according to equation l|86j l. The bold line is the theoretical prediction. 
Numbers in the graph denote avalanche duration in ms. (Adapted from [113]). (Reprinted with 
permission from [113]). [113] Copyright 2004 Elsevier. 




0.1 1 10 

duration 7~(ms) 



Figure 26. Average avalanche size as a function of its duration for polycrystalline FeSi material. 
Data are from the same experiment of figure I25I The solid line is the predicted scaling S ~ T 2 . 
(Reprinted with permission from [5]). [5] Copyright 2005 by Academic Press. 

The signal-signal time correlation function is denned as 

C(t) = ((v(t)-c/k)(v(0)-c/k)). (89) 

By deriving the probability R with respect to t and using the Fokker Planck 
equation, it is easy to get a differential equation for C (see [31] for the details), 
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which turns out to be 

dC{t) 



dt 



kC(t) . (90) 



Therefore C(t) = c/fcexp (—kt), which corresponds to a Lorentzian power 
spectrum 

(91) 

Equation (|9ip captures the correct u>~ 2 large frequency decay found in exper- 
iments, however it saturates to a constant for small u, without showing any 
peak. 

To get a better prediction of the power spectrum, let us go back once again 
to the original formulation of the ABBM model, which includes a more general 
case, that we have not considered up to now, in which a finite correlation length 
£p is associated to the pinning field [31]. A finite £p has the effect of confining 
the pinning field, and physically measures the length over which the pinning 
field decor relates. 

In most cases it can safely be assumed tip — * oo, as done in most literature, 
and in the analysis reported in the previous sections. However, the existence of 
a finite correlation length in the pinning field affects the low frequency behavior 
of the power spectrum. The power spectrum for the case of a finite pinning 
correlation length is calculated in [31] in a high field velocity approximation, 
and results to be 



{oJ 2 + k 2 )(uj 2 + T p 2 ) 



F (") = , o . . -2, > ( 92 ) 



where rp = £p/c is a characteristic time associated to the correlation length 

Equation ([92]) reproduces the low frequency part of the spectrum. It predicts 
the existence of a peak at lom = (kTp) 1 / 2 ~ c 1 / 2 , with a maximum value of 
the spectrum Fm = , -i^ which scales as c for small driving rates, as 

observed in experiments. 

It would be interesting to understand how the presence of a finite tip affects 
other properties of the avalanches, as the decorrelation of the pinning field at 
some finite length scale might be responsible of some of the observed deviations 
from the tip — > oo ABBM theory in large events. 

The deviation from the 1/ui 2 prediction for the tail, which is observed at 
intermediate frequencies in some experiments might have the same origin as 
the deviation from scaling observed in the pulse shape scaling in the previous 



February 18, 2009 15:51 



Advances in Physics adpl.21 



Barkhausen noise: mean field and beyond 61 




Figure 27. Comparison of the power spectrum F(uj) with the average avalanche size (S) as a 
function of the inverse duration for a polycrystalline FeSi material. The solid line is the predicted 
scaling S ~ T 2 . (Adapted from [117]). (Reprinted with permission from [117]). [117] Copyright 

2002 by Elsevier. 

subsection [117]. This is suggested by the comparison of the power spectra with 
the average size versus inverse avalanche duration. The two curves deviate for 
the same value of u (see figure |27"|) . The origin of the deviation in the size 
versus duration plot observed in some experiments is, however, still unclear. 

A possibility is that the eddy current retardation that we will discuss in the 
next section is responsible for the intermediate frequencies deviations from the 
asymptotic scaling, as suggested in [118] based on simulation of the modified 
ABBM equation. We will come back to this in the next section. 

4.9 Summary of ABBM predictions and comparison with the 
phenomenological observations 

Let us summarize the ABBM predictions and compare them with the phe- 
nomenological observations from BN experiments on soft materials with long 
range interactions: 

(i) The distribution of the amplitude of the BN signal follows a power law 
decay with an exponent 1 — c as observed in experiments (see equation [5]) . 

(ii) The distributions of sizes and durations of avalanches can be calculated 
in terms of recurrence properties of the process described by equation [MJ 
and turn out to be power laws, with exponents a = 2 and r = 3/2 in the 
quasi-static limit c = 0, as found in experiments. 
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(iii) Equation [M] predicts a linear dependence of the exponents a and r on the 
driving rate c, given by a = 2 — c and r = (3 — c)/2, the same obtained in 
experiments. 

(iv) There exist a threshold c = 1 above which the process in equation [M] 
become transient. This corresponds to the threshold in the driving field 
rate above which the BN pulses collapse into a unique large avalanche of 
magnetization reversal. 

(v) The bias term in equation [34"l deriving from the demagnetizing field, is 
responsible for the exponential cut-offs in the distributions of sizes and 
durations of avalanches. The dependence on k of the cut-off predicted by 
the ABBM equation is however different from that observed in experi- 
ments. 

(vi) The average excursion of a trajectory for the stochastic process in equation 
1341 can be calculated and is related to the average pulse shape in BN. 
However, the predicted shape is symmetric, while experimental pulses show 
a clear leftward asymmetry. 

(vii) The power spectrum decays at large frequencies as u~ 2 . A finite correla- 
tion length for the pinning field produces a peak at small frequencies at 
•jJm ~ c 1 / 2 of amplitude Fm ~ These results are compatible with experi- 
mental ones, however, deviations from the predicted behavior are observed 
at intermediate frequencies in some experiments. 

These predictions can be compared point to point with the experimental 
observations summarized in sec 12.4.91 Most of the phenomenology is captured 
by the ABBM model, however, two main unexplained discrepancies between 
experiments and theory appear: one is the scaling of the cut-offs with the 
demagnetizing factor k, the other is the presence of a leftward asymmetry in 
the pulse shape. 

The first inconsistency can be dealt with by going back to the full elastic in- 
terface model, as we will discuss in briefly the next subsection. The second one 
is more serious, as the pulses asymmetry is related to non-universal properties 
of the magnetization reversal of microscopic origin, and we will postpone its 
discussion to the next section. 



4.10 Cut-off out of mean field approximation 

Let us go back for a moment to the elastic interface model discussed in section 
13.41 This model has been extensively simulated both in the infinite range case, 
corresponding to the mean field, and in finite dimension d. The case d = 3 
is supposed to give results that, apart from logarithmic corrections, coincide 
with the mean field, since d = 3 is the upper critical dimension for this system. 
This is confirmed in simulations, where one gets in both cases, for distribution 



February 18, 2009 15:51 Advances in Physics adpl.21 



Barkhausen noise: mean field and beyond 63 

of sizes and durations of the pulses, the exponents a and r predicted by the 
ABBM theory, both in the adiabatic limit c — > 0, and for finite driving c. 

However, the scaling with k of the cut-offs Sq and To in the distributions is 
different in the case d = 3 and in mean field. The simulations on the infinite 
range model confirm the results of section 14. 6t they predict that the cut-offs 
scale as Sq ~ k~~ 2 , and To ~ k , consistently with the results obtained for the 
ABBM model. However simulations in d = 3 are not compatible with these 
results. 

Earlier simulations [56] in d = 3 give Sq ~ k~ 1 , and To ~ A; -1 / 2 , both in the 
quasi-static limit and at finite c. In [56] the full size and duration distributions 
are fitted by 

P(S,k)^S-^ 2 fs(S/k-^), (93) 

and 

P(T,k) ~ T- ( - 2 -^f T (T/k-^ T ) , (94) 

with 7s ~ 1, and jt — 1/2, which are claimed to be compatible with the 
scaling obtained in the experiments reported in the same paper, where the 
cut-offs dependence on k has been investigated by performing experiments on 
samples with different geometries. It has to be taken into account the extrap- 
olation of a precise cut-off value is complicated by the further dependence of 
the distributions on the driving field rate c. 

In more recent papers, different values are reported [54,57], corresponding 
to 7s ~ 2/3, and jt — 1/3. These latter results are also supported by the 
following argument. The maximum interface jump is limited by a correlation 
length £, which is the length at which the dipolar interaction term, of order 
overcome the restoring force. In d = 3 (where the interface is not rough), 
this is of order /c£ rf_1 , where ^ d ~ 1 h is the size of the avalanche. This gives 
£ ~ k~ x l d , and therefore S ~ I 4 - 1 ~ Ar^ 1 )^ = A; 2 / 3 . 

In the same papers [54,57], results from more recent experiments are also 
reported, corresponding to 7s ~ 0.6, and 7t ~ 0.3, that could be compatible 
with the values 7s — 2/3, and 7t ~ 1/3. In these works, the experimental esti- 
mates of the dependence on the demagnetizing field is extracted by repeating 
experiments on the same sample, which is progressively cut in order to change 
the value of k. This procedure allows to change the value of k while stress and 
internal disorder are kept constant in the repeated experiments (see figure [T3| 
which also reports the results of experiments on materials belonging to the 
short range universality class not discussed here). The experimental results 
in [54,57] are therefore expected to be more reliable than previous ones, where 
different samples where used as in [56]. 
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The scaling with k of the cut-offs in the distributions of sizes and durations 
of the avalanches seems to agree with that predicted by the elastic interface 
model at the upper critical dimension d = 3, although this prediction is dif- 
ferent from the mean field one. The origin of the discrepancy between mean 
field behavior and the behavior at the upper critical dimension is not fully 
understood, and we believe that this point would still need further analysis in 
order to be definitely clarified. 



5 Beyond ABBM: Asymmetry in the pulse shape 

As we have seen in the previous section, despite its extreme simplicity, the 
ABBM model is able to capture a number of phenomenological evidences from 
BN experiments. However, the analysis on the scaling functions is more sub- 
tle: although the average pulses do approximately collapse onto each other 
once rescaled by their duration, by comparing the average shape of rescaled 
experimental Barkhausen pulses with the prediction from the ABBM model, 
one clearly observes a discrepancy: while the ABBM model predicts a sym- 
metric shape, experimental pulses are systematically characterized by a visible 
leftward asymmetry, as can be seen in figure [241 The leftward asymmetry in- 
dicates that, unlike what would be expected from standard inertia, avalanches 
start fast and end slowly. 

Given the success of ABBM theory in reproducing many other aspects of 
phenomenology, this inconsistency has been considered quite enigmatic, and 
the issue of understanding its origin has remained an open problem for some 
time. 

The first attempts to identify the origin of the asymmetry have been done 
in [114], and more extensively in [115]. In these works, the authors pose the gen- 
eral problem of analyzing the behavior of the average excursion in a stochastic 
process, in order to understand to which extent this quantity is universal, and 
what kind of information it gives on the process. The idea was that this anal- 
ysis could suggest in which direction the ABBM equation had to be modified 
in order to recover a leftward asymmetric pulse shape. 

For a large class of stochastic processes, a scaling law of the form [7T] is 
obeyed, in an appropriate time regime, with a scaling function that is, to a large 
extent, independent of the details of the single increment distribution, while 
it encodes interesting information on the presence and form of correlations. 
Indeed the introduction of correlated noise, is found to induce asymmetric 
shapes. However, in the cases considered in [114,115], the scaling functions are 
characterized by a negative skewness, corresponding to a rightward asymmetric 
shape. 

More recently [118, 119], the origin of the asymmetry has finally been identi- 



February 18, 2009 



15:51 



Advances in Physics 



adpl.21 



Barkhausen noise: mean field and beyond 



65 



fied as a transient effect of eddy currents. Eddy currents are generated in con- 
ducting magnets as a response to domain wall displacements. This response, 
however, is not instantaneous and acts as an anti-inertial effect on the do- 
main wall dynamics. It turns out that the effect of this retardation can indeed 
be accounted for, to a first approximation, by including a negative effective 
inertial term in the stochastic equation of motion for the domain wall. 

5.1 Experimental evidence of asymmetric pulses 

Figure shows an attempt to collapse average pulses of different durations, 
obtained from a BN experiment on a polycrystalline FeSi sample, onto the 
same function. The pulses do approximately collapse, when properly rescaled, 
however, the supposedly universal shape is clearly not symmetric, as would be 
predicted by the ABBM model (see figure l24l) . The asymmetry appears to be 
more pronounced for small avalanches. 

To get a quantitative measure of the asymmetry the skewness can be calcu- 
lated as: 



where t = dt{v(t,T))/T. 

Figure [25] shows the skewness as a function of the pulse duration in an 
experiment on a partially crystallized Fe^Co2iBi^ ribbon. The skewness is 
always positive, corresponding to the leftward asymmetry. It shows a peak for 
Tp ~ 400/US, and decays to zero as T goes to infinity. 

The fact that, except for the initial increase at very small durations, the 
asymmetry becomes smaller and smaller for larger and larger avalanches sug- 
gests that the ingredient still missing in the ABBM theory must involve mi- 
croscopic details of the system that become negligible on larger scales. The 
presence of a peak indicates the existence of some microscopic timescale. 

5.2 Retarded eddy current pressure 

To understand why ABBM is unable to reproduce the leftward asymmetry of 
the average pulse shape, we have to go back and question the assumptions 
which were made in its derivation. The equation of motion for the magnetic 
interface was derived by assuming that the work done by the total effective 
field (the applied field corrected by the demagnetizing and pinning fields) 
balances at any time the energy dissipated by eddy currents, which, in the 
quasistatic approximation, is proportional to the instantaneous velocity. The 
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Figure 28. 



Skewness as a function of the avalanche duration measured in an experiment on a 
partially crystallized Fee4Co2iBi5 sample. (Reprinted from [119]). 



equation is overdamped, and neglects inertial effects. Let us mention that the 
existence of inertial effects in domain wall dynamics, has been demonstrated 
by Doering [120], who showed that a moving wall differed in energy from the 
wall at rest by a term due to gyromagnetic effects and proportional to the 
square wall speed. He identified this as a kinetic energy term and defined the 
constant of proportionality as one-half of the mass of the wall. The so called 
Doering mass Mjj results to be proportional to the inverse wall width divided 
by the squared gyromagnetic ratio (for a simplified treatment see [121], and 
also [122] for a review). Inertial effects are therefore always present in domain 
wall dynamics, however, they are not relevant up to the GHz band, and can 
be safely neglected in the equation of motion for the magnetic wall in other 
regimes, where eddy currents damping largely overwhelms the contribution 
from Doering mass. 

The assumption that the work done by the net field is at any time propor- 
tional to the instantaneous velocity is, however, only approximately correct: 
eddy currents take a finite time to set up after the corresponding wall dis- 
placement, and also, they persist for a finite time after the magnetic reversal. 
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Because of this delay, the eddy pressure on the wall at a given time t of the 
avalanche depends on the history of the motion: is not strictly proportional 
to the instantaneous velocity, but rather it is a weighted average of all the 
velocities from the beginning of the avalanche up to time t. 

To go beyond the quasistatic approximation, one has to investigate the dy- 
namical effects eddy current dissipation. Following the lines indicated by [123] 
this has been done by calculating the magnetic field generated by the eddy 
currents from the full Maxwell equations [119,124]. 

Consider a sample with dimension x 6 [—a/2, a/2], y G [—6/2,6/2], and 
infinite in the z direction, divided in two magnetic domains by a rigid domain 
wall on the yz plane, that starts to move from the zero-magnetization position 
x = 0, as in figure 1291 

The displacement of the magnetic wall in the conducting medium induces 
a flow of eddy currents that generates a magnetic field H e , which, in this 
geometry, is parallel to the z axis: 



Neglecting the displacement currents with respect to the ohmic currents, the 
Maxwell equation for H e is 



where a and \i are the electric conductivity and the magnetic permeability of 
the medium respectively. 

Equation ((971) is a diffusion equation for the eddy field, with a typical 
timescale for diffusion proportional to o[i. In the quasi-static approximation 
the permeability is assumed to be negligible within the magnetic domains, and 
equation ((971) reduces to V 2 H e = 0, corresponding to an instantaneous eddy 
field propagation. 

This approximation is correct as long as there is a clear separation of 
timescales, which happens when the observed quantities do not vary apprecia- 
bly in time on a scale of the order of the characteristic time for diffusion. In 
that case, one can assume an instantaneous response of the eddy field as given 
by the static Maxwell equation. However, this is not the case for Barkhausen 
avalanches, especially for short ones, and therefore the dynamic effect of eddy 
field has to be taken into account. 

Let us go back to equation (f97j) . which has to be solved subject to the 
appropriate boundary condition 



H e = H e (x,y,t)z. 



(96) 




(97) 



H e = 



(98) 
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Figure 29. Horizontal arrows indicate the directions of the magnetization in the two domains. The 
black vertical arrow indicates the direction of motion of the wall. 

on the sample surface, and with the Faraday condition 

d x H e (0 + ,y,t) - H e d x H e (0- t y,t) = 2aM s v(t) , (99) 

that regulates the discontinuity of the field across the wall. In equation [99] M s 
is the saturation magnetization, and v(t) is the instantaneous velocity of the 
domain wall. 

The solution of equation (|97p has been derived in [124] , by expanding H e in 
its Fourier components 



H e (x,y,t) 



1 



dujF(x, y, u)e 



iujt 



f 2ir j-oo 
Each component has to satisfy 

VF(x, y, u) = r 2 F(x, y, u) , 



(100) 



(101) 



with r 2 = iuiafi. The condition for the magnetic eddy field to be zero on 
the sample boundary implies that F(x,±b/2,uj) = F(±a/2,y,Lu) = 0. The 
solution of equation (|118p for the Fourier components with these boundary 
condition is of the form 



F(x, y,uj) = A n(x, u) cos(A n y) 



(102) 



n=0 



where A n satisfies 



d%A n (x,uj) = AlA n (x,Lj) 



(103) 
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with A 2 = A 2 + r 2 , to be solved separately for x > and x < 0. The condition 
on the sample boundary in the y direction implies cos(±6/2A ra ) = which 
fixes X n = (2n + 1)tt jb. The condition on the other boundary in satisfied by 
choosing 

A n {x, u) = C n {uj) sinh(A n (| x | -a/2)) , (104) 

so that 

CO 

F(x,y,u) = ^C n (w)sinh(A n (o;)(| x \ -a/2)) cos(A n y) . (105) 

n=0 

The functions C n (oj) are fixed by the Faraday condition around the wall, 
which in Fourier space reads 

d x F{0 + ,y,uj) =-d x F(0~,y,uj) = aM s v(u) , (106) 

where v is the Fourier transform of the velocity of the wall v(t). Equation 
(1106j) implies 

CO 

C n {oj)A n (uj) cosh(A„(u;)a/2) cos(A n y) = aM s v(cv) . (107) 

n=0 

Multiplying by cos(A m y), integrating in [—6/2, 6/2], and using the orthogonal- 
ity relations f^{/ 2 dy cos(A n y) cos(A m y) = S n>m b/2, and f^ 2 dy cos(A m y) = 
(-l) m 2/A m one gets 

C n (u) = (-ir ia f s 1 — v(u) , (108) 

b \ n k n [uj) cosh(A n (u;)a/2) 

so that the solution for the Fourier components of the field H e is finally found 
to be given by 

„, . 4aM s A, ,_ sinh(A n (u;)(| x I — a/2)) » . . 

F(x,y,u) = — -Jl^-i" a 7 \ u\ ( / /9 > ^ • 109 
b A ri A n (cj)cosh(A n (cj)a/2) 



To calculate the pressure on the wall due to the presence of eddy currents, 
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only the value of the field at x = is needed, where F simplifies to: 



cvn \ W ^n ^Mi tanh(A n (u;)g/2) 
F(0,y,w) = 2^(-l) — t A 7-n cos(\ n y)v(uj) . (110) 

n=0 n n ^ ' 

The average eddy current pressure on the wall is then obtained by integrating 
the magnetic field over y at the wall position x = 0: 

21 f h l 2 21 1 f°° f h l 2 

Pit) = — / dyH e (p,y,t) = — -= / due"* / dyF(0,y,u). 

b J-b/2 b V2VT J-oo J-6/2 

(111) 

In real space the pressure at time t is given by a convolution of velocities of 
the wall at all times prior to t with some response function /: 

P(t) = -L f° d«7(t - *)/(*), . (112) 

V27T J-oo 

which, in Fourier space, corresponds to 

P{oj) = -v(u)f(u) . (113) 

Replacing F(0,y,u) in equation (jllip with its expression (jllOp . and perform- 
ing the integral over y, we get, for the Fourier transform of the response func- 
tion, the expression: 

IfiPa ^ tanh(A,,Ha/2) 

/m =— E • (m 

To proceed further it is convenient to simplify the u dependence in equation 
(|114[) . which can be done by using the relation 

OO j 

tanh(z) = 8z g ^ {2k + l?+Az 2 > ( 115 ) 

with z = A n (u;)a/2, which gives 

* _ 64al 2 cr ~ l 

/M -^i A 2 (^ + a^H)- (116) 
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Replacing A n and A n with their expressions, one gets 



QAI 2 a 
ab 2 e 2 



oo 



1 



E 



n 



2 uj (k 2 to a + n 2 u>b + iui) 



(117) 



n,fc=l 



where Yl indicates a summation over odd numbers only, and u a = r^ 1 = 
Tr 2 /ea 2 , uj b = t^ 1 = tt 2 jeb 2 . 

From equation (|117|) . the expression for the response / in real space is ob- 
tained by anti-transforming term to term in the double sum, and replac- 
ing each term of the form \/{ujq + ioj) with its inverse Fourier transform 
\/2vr exp(— u>ot)8(t) (where 9(t) is the Heaviside theta function). This leads 



where 0Jk,n = T kn = k 2 uj a + n 2 ujb. The form (|118p of the response function gives 
through the convolution (lllip the explicit expression of the retarded pressure 
on the wall. 

The response function finally results to be the sum of simple exponential 
relaxations, with different relaxation times. The most relevant relaxation time 
is the largest one, which is 



This is the quantity that has to be compared with the typical timescale of the 
phenomena under study to decide whether or not the dynamical effects due 
to eddy current relaxation are negligible. 

In a typical Barkhausen noise experiment, t^o is of the order of 5/is. Al- 
though this time is small with respect to overall duration of a Barkhausen 
pulse, the magnetization within an avalanche does vary appreciably on that 
timescale (see figure [H]). 

5.3 Generalized ABBM: negative effective mass 



to: 




(118) 




(119) 



In terms of the equation of motion, the effect due to the non-instantaneous 
response of the eddy field to the wall displacement is accounted for by replacing 
the left hand side of equation (129j) with the retarded pressure P(t). This leads 
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to a modified ABBM equation: 



t 

dsf(t- s)v(s) = H{t) - km + W(m) , (120) 

o 

that includes the time nonlocal effects. 

The eddy dissipation keeps memory of the past velocities of the interface, 
within a timescale of the order of the eddy field diffusion relaxation time. 

To clarify the effects of the response function, it is useful to consider the case 
of a simple exponential decay f(t) = T/tq exp(— t/ro), which corresponds to 
the contribution of a single term in equation (I118p to the equation of motion. 
As long as ro is small with respect of the duration of the avalanches considered, 
only velocities at small times before t will contribute to the convolution, so 
that we can expand v (s) around t and perform the integrals: 

ft e~ s l T ° 

T / ds v(t - s) ~ Tv(t) - TT v'(t) . (121) 

Jo T o 

A simple exponential response function would therefore modify the equation 
of motion by introducing an inertial term 

Tv(t) + Mv'(t) = H(t) - km + W(m) , (122) 

and would lead to the identification of an effective mass M = — IYo, which 
turns out to be negative. 

When the full response function is considered, a series of mass and damping 
terms are generated, each coming from a single exponential. Taking only the 
leading order in /u, this procedure leads in a similar way to the identification 
of a damping coefficient 

a7r 4 

and of a negative effective mass 

M = - V Mb/a), (124) 
an 

where Si (a) = J2 n ,k iP n'+a'fc" and S 2(a) = £„ >fc ^ (n2+ ^ fc2)2 , and £ indi- 
cates a summation over odd numbers only. 
As in the case of the simple exponential relaxation, the ratio between mass 
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and damping defines a characteristic time 



t = \M\/T 



= u, 



,-i Sa(V°) 
b Si(6/a)" 



(125) 



Keeping all orders in [i results in a frequency dependent effective mass, that 
has been calculated in [124] and is given by 



The effective mass turns out to be negative at all frequencies. 

The leftward pulse asymmetry observed in Barkhausen experiments is indeed 
consistent with a negative effective mass: the avalanches start fast and end 
slowly, which is exactly the opposite of what one would expect from standard 
inertia. Indeed, since the retarded pressure at time t is given by a weighted 
average of previous velocities of the wall up to time t, this anti-inertial effect 
can be understood by observing that the effective average velocity is smaller 
than the instantaneous one for an accelerating wall, while the opposite is true 
when the wall decelerates. 

5.4 Effective mass dependence on sample geometry 

The role of sample geometry on eddy current retardation has been studied 
in detail in Ref. [124]. It turns out that the only geometrical parameter that 
significantly affects the damping coefficient, effective mass, and characteristic 
time, is the smallest sample dimension. 

Figure l30l shows the variation of damping coefficient T, characteristic time 
r, and of the modulus of the effective mass M, as a function of the aspect 
ratio r defined as 



and keeping constant the section area A = ab. These three quantities have a 
maximum in the case of a square rod. Therefore this is the geometry where, 
A being constant, the maximum asymmetry is expected to be observed. 

Figure EU shows the variation of the the same three quantities as a function 
of the most relevant parameter, which is the smallest sample dimension d = 
min(a, b). The damping coefficient T, characteristic time r, and modulus of 
the effective mass M, increase as d, d 2 , and d 3 respectively, and saturate to 
the values A 1 / 2 , A, and A 3 / 2 , respectively. The role of the two dimensions a, 




(126) 



r = max(a, b)/ min(a, b) 



(127) 
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Figure 30. Damping coefficient T, characteristic time r, and modulus of the effective mass M, as 

functions of the aspect ratio r. 



and b, transverse and parallel to the wall is however not symmetric, and the 
effect of eddy current retardation is expected to be more severe, other things 
being equal, when the transverse dimension is the smallest. 

The analytical expressions of T, r, and M used in figures [30l and [31] are 
given in [124]. 

It would be interesting to repeat BN experiments in changing sample ge- 
ometries in order to test these results. 

Let us comment on the relevance of the eddy inertial term when compared to 
the usual inertial term which is associated to the Doring mass Md . The values 
of Md are about four orders of magnitude smaller then the eddy mass M in 
typical sample geometries used in BN experiments. Indeed the Doring mass 
originates from a gyroscopic precessional effect, whose associated time scales 
are much faster then those relative to the eddy mass. However, the analysis 
reported in [124] on the dependence on sample geometry predicts a fast decay 
<i~ 3 of the mass with the sample thickness, which means that the two masses 
may become comparable in thick films. 
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5.5 Comparison between theory and experiments 

By simulating numerically the modified ABBM equation (|120p one observes 
that, whereas the distributions of avalanche duration and size are unaffected 
by the addition of the inertial term, the pulse shapes become asymmetric and 
bear a remarkable similarity with the experimental ones (see figures l32| [33| 
where the curves are normalized with N = / dt(v(t, T))/T) [119]. Note that 
eddy currents only change the internal dynamics of the avalanches, so that 
at least the distribution of avalanches sizes' is expected to be left unaltered, 
while correction could be present in the power spectrum (see below) and in 
the distribution of avalanche durations' (not observed). 

The skewness of the simulated model results positive, indicating a leftward 
asymmetry, corresponding to a negative mass, while it would be negative, in- 
dicating rightward asymmetry, in the standard case of a positive mass. The 
qualitative behavior of the skewness as a function of the avalanche duration 
in correspondence to a negative mass (see figure is similar to the one ob- 
served in the experiments (see figure l28j) : it is always positive, it increases from 
zero, shows a peak at some characteristic time, and it decays to zero for long 
avalanches. The existence of a peak can be used to extract the characteristic 
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Figure 32. Asymmetry of the rescaled pulse shape: pulses obtained from Barkhausen noise 
experiments in a partially crystallized Fe%^Co2\B\^ ribbon. The shapes for different durations T. 

(Reprinted from [119]). 

relaxation timescale, which corresponds to the ratio between mass and damp- 
ing constant, and to make a quantitative comparison between the model and 
experiments. 

A further quantitative test of the theory against experiments includes the 
comparison of the location of the peak Tp that characterizes the skewness ver- 
sus duration plot. Skewness extracted from experimental pulses shows a T p at 
a value which is roughly twice that obtained from simulations of the modified 
ABBM equation with the microscopic parameters of the corresponding exper- 
iment. This may be due to the fact that the model considers just one single 
domain wall, while many domains may be present in general in a sample. It is 
indeed expected [123] that the relaxation time increases with the number of 
domains. 

Another interesting point to underline is the effect of the negative inertial 
term on the power spectra. Power spectra obtained by simulating the modified 
ABBM equation show a deviation from the large frequency asymptotic regime 
at a frequency uj* = Tp 1 very similar to the one observed in some experiments. 
Moreover, experimental data also show a correspondence between the peak in 
the skewness and the frequency at which the deviation appears in the spectrum 
[118]. 

Summarizing the discussion of this section, we can conclude that the left- 
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Figure 33. The shapes obtained from the modified ABBM model with the parameters 
corresponding to those of the experiment in figure l32l (Reprinted from [119]). 

ward asymmetry of pulse profiles in BN is caused by the non instantaneous 
response of the eddy field to the wall displacement, that can be accounted for 
by associating a negative effective mass to the moving wall. The asymmetry 
is a non-universal property, that depends on the material parameters a, /i, 
and M s and on the geometrical dimensions a and b of the sample, and there- 
fore it allows to extract information about the characteristic time ruling the 
microscopic dynamics. 

Since the asymmetry is related to the existence of this microscopic timescale, 
its magnitude depends on the avalanche duration: for avalanches long enough 
with respect to the characteristic time, the microscopic details of the mate- 
rial and the geometry of the sample become irrelevant, pulse-shapes become 
symmetric, and universality is expected to hold. 

Let us finally comment on the physical meaning of the eddy inertial term. 
It has always to be kept in mind that the negative effective mass is just a con- 
venient way to account, to the first order, for the history dependent damping. 
The effective mass, however, is unable to describe for example the very short 
time dynamics, when the domain wall starts to move from its position at rest. 
Indeed, in this case, the negative mass would cause the wall to move backward, 
which is certainly not the case. Here the results of this section should be used 
with more care: when the wall starts to move there is no previous motion his- 
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Figure 34. Skcwncss from simulations of the modified ABBM model. Filled symbols correspond to 
positive mass (negative skewness), empty symbols to negative mass (positive skcwncss). (Reprinted 

from [119]). 

tory, and therefore both the eddy damping and the eddy inertial term are zero. 
In the very short time dynamics the Doring mass, which is positive, represents 
the most relevant inertial term in the equation of motion. Therefore standard 
inertia governs the dynamics, until damping and eddy mass build up. 



6 Other crackling systems: granular materials, seismic activity in 
earthquakes 

Crackling noise emerges in very different contexts, and the similarity in how 
completely different systems behave suggests the existence of some common 
principle in the underlying physics. As an example of how ideas developed 
in the study of one crackling system can stimulate advances in other fields, 
we report here some analogies between BN and two other crackling system: 
sheared granular materials and seismic activity in earthquakes. 
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Figure 35. Comparison of power spectra where deviation from 1/uP is observed at intermediate frequencies with the one from simulations, (a) 
Power spectrum from a simulation of the modified ABBM model. (Reprinted with permission from [118]). [118] Copyright 2007 by Elsevier, (b) 
Power spectrum for a partially crystallized Fe^4Co2iBl5. (Reprinted with permission from [81]). [81] Copyright 1996 by Elsevier, (c) Power 
spectrum for a polycrystalline 7.8% FeSi. (Reprinted with permission from [58]). [58] Copyright 1997 by World Scientific, Singapore. 
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6.1 Dynamics of granular materials under shear 

A granular medium, subject to a slow loading rate, has an intermittent shear 
response characterized by large fluctuations, typical of crackling systems [125, 
126]. The large fluctuations in this case originate from the rearrangement of 
the network of force chains along which the stress propagate. 

Most of the past theoretical activity on granular matter has focused on 
average properties, that could allow the formulation of macroscopic laws gov- 
erning the dynamics. However, more recently it became clear that a deeper 
understanding of the mechanisms that regulate the dynamics on smaller scales 
and cause the large fluctuations, cannot be reached without going through a 
detailed analysis of the slip statistics. 

In this framework, a very interesting model has recently been proposed to 
interpret a series of experiments on granular materials sheared in a Couette 
geometry [7]. This model, is based on an assumption similar to the one that 
underlies the ABBM model for BN, namely that the resultant from the forces 
acting on the system performs a Brownian motion. 

In the experiment a cylindrical cell is filled with monodisperse glass beads 
with 2mm of diameter. The shear is created by rotating an annular plate over 
the top of the circular channel, driven by a torsion spring with stiffness k. 
After being initialized to reach the stationary state, the system is then run 
at a slow driving velocity v. The experimental device allows to measure the 
angular position of the plate, and the deflection of the torsion spring. In a 
given range of driving velocities, the system displays a stick-slip regime, where 
repeatedly the torque on the plate accumulates as the spring winds up, until 
the plate slips. 

The instantaneous velocity of the plate 6 as a function of time, shows the 
typical intermittent behavior of crackling signals, with pulses of largely fluctu- 
ating amplitudes. A statistical analysis of the signal similar to the one usually 
done for BN emissions reveals the presence of power law distributions of the 
slip events in the stick-slip phase, as it was already observed in a previous 
experiment [126]. 

In this case the pulse size s corresponds to the angular variation during the 
slip event. The sizes and durations distributions show a more complex shape 
than in the BN case, since the initial power law decay is followed by a peak 
at larger scales, in correspondence to some cut-off value. The position of the 
peak is independent of the driving velocity v. 

As in the case of BN, the intermittent stick-slip phase is observed only as 
long as the driving velocity is not too large. Above a threshold value, a tran- 
sition to a steady sliding regime occurs, that is analogous to the transition 
between the intermittent BN noise regime, and the unique avalanche magne- 
tization reversal that occurs at high sweep field rates. 
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In ref. [7] the authors propose a simple model that quantitatively reproduces 
the phenomenology, and which is based on some minimal hypothesis. The 
equation of motion for the stick-slipping plate is obtained by equating the 
inertial term to the driving force on the spring diminished by the counteracting 
friction of the granular media. 

I0 = k(vt-9)-F. (128) 

Here I is the moment of inertia of the system, and the whole complexity of the 
granular matter dynamics that give rise to the wide fluctuations is accounted 
for in the frictional torque term F. 

Baldassarri et al. assume that, as a first approximation, the instantaneous 
frictional torque is a function of 9 and 9 only (since F is constant when the 
disk is stuck, there is no direct dependence on time t in F), and that F can 
be split into a deterministic part Fd{6) plus a random contribution F r {6). 

The deterministic part of the instantaneous friction F^O) is obtained by 
fitting the experimental average velocity dependence of the frictional torque 
with the expression 

F d {9) = F + 1 {9- 2v m(l + 9/v j) , (129) 

where Fq is the average static friction force, vq is the minimum average torque 
and 7 is the high velocity damping. Note, that the results obtained, however, 
do not depend very much on the precise choice of the function Fj. 

The fluctuating part of the instantaneous friction is the one that accounts for 
the continuously rearranging disordered structure of the force chain network 
present in the granular medium. This is assumed to be described by a confined 
Brownian process 

^l = v (9)-aF r (9), (130) 

where r\ is a Gaussian quenched uncorrelated noise, and a is an inverse correla- 
tion length. The idea is that the frictional torque changes of a random amount 
after every small slip, but eventually decorrelate on large slips involving a com- 
plete rearrangement of the grains, resulting in a confinement of the Brownian 
process. The quantity a -1 plays the same role as the correlation length of the 
pinning field £p introduced in the ABBM model (see section . 

This theory shows a number of analogies with the ABBM theory for BN. 
The granular matter dynamics is more complicated than the one of magnetic 
domains, due to the presence of inertial effects and of a shear rate weakening 
torque. However, both theories are based on the same assumption that the 
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effective force acting on the system performs a Brownian motion. Indeed, the 
presence of Brownian correlations in the effective force evolution, due to col- 
lective effects (the sum over many pinning forces over the domain wall in one 
case, the formation and destruction of force chains in the granular media in 
the other) may be of more general validity in driven dissipative systems. 

The equation of motion is also similar in the two cases, and, although, 
the one describing the granular material includes more complex effects, it 
reduces to the ABBM equation in the limit /, vq — > of vanishing inertia 
and minimal average torque. In this case, or in an appropriate velocity regime 
vq *C v <C 7//, where the additional effects can be neglected, all the results 
obtained for the ABBM model are valid, once magnetization is interpreted 
as a slip, the magnetic field rate as a shear rate, the demagnetizing factor as 
the stiffness of the torsional spring, the eddy current damping as a frictional 
damping, and the random pinning as a random friction. 

The analogy also suggests that the transition from the stick-slip phase at 
slow drivings to the sliding phase at fast drivings might possibly be explained 
also in this case as a transition from transient to recurrent for the correspond- 
ing stochastic process. 



6.2 Seismic activity in earthquakes 

Another kind of crackling noise that shows analogies with BN is the seismic 
activity during earthquakes. 

Earthquake phenomenology exhibits a number of power law distributions 
including the Gutenberg-Richter frequency-size statistics [127, 128] and the 
Omori law for aftershock decay rates [129,130]. 

In literature, the area of the region of the fault that slips in an earthquake is 
called the moment. The probability of occurrence of an earthquake of a given 
moment M decays as a power law according to 



The exponent f3 is the same, within errors, for most sets of data, collected in 
different times and regions. 

As in the case of BN, scaling functions have been introduced as additional 
interesting observables, which provides a powerful tool both to check univer- 
sality, and to test the theory against data. In particular, the average moment 
rate profiles at fixed total earthquake duration, have been studied in [30], and 
compared with the prediction of a mean field theory for earthquakes dynam- 
ics [131]. 



P{M) ~ Af-( 1- fl . 



(131) 
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Figure 36. A collapse of averaged earthquake pulses' shapes, with a duration of T, with two to ten 
earthquakes averaged for each value of T. The collapse was obtained using the mean field theory 

for which the universal scaling function is claimed to be g(x) = Ax(l — x) with x = t/T [51]. Note 
the apparent asymmetry to the left in the observed data while the theoretical curve is symmetric 
around its maximum. In the inset the raw data and the averaged data before collapse are shown. 
(Reprinted with permission from [30]). [30] Copyright 2006 by the American Physical Society. 

The average moment rate is supposed to behave as 

(dM(t, T) /dt) ~ Tg(t/T) . (132) 

In paper [30] the universal function g(t/T) is claimed to be: 

g(x) = A(x(l - x)) . (133) 

This result was obtained by mapping the model onto the mean field model for 
BN, for which the result of g being an inverted parabola was claimed in [51]. 

It is clear from figure [36] that in seismic data the moment rate of a given 
duration T increases faster then it drops, as the corresponding average moment 
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shapes show a leftward asymmetry similar to that observed in BN, while, also 
in this case, models predict a symmetric shape. 

Seismic movements are caused by the motion of fault planes in response 
to the accumulated stress. The analogy with BN suggests the possibility that 
the origin of the asymmetry is similar in the two cases. The asymmetry may 
indeed originate from some non-local dynamical effect, that, in the case of 
earthquakes could be caused by the presence of stress overshoots [132]. It 
would be interesting to see how the skewness of the average pulses changes 
with the duration of the seismic event (which is not clear for the data in 
ref. [71]), as this can help to understand whether or not also in this case the 
asymmetry disappears in large events, and whether universality is recovered 
oon large time scales. 

In the case of earthquakes, the analogy with domain wall dynamics is even 
more stringent than in the case of sheared granular media. Indeed, the mean 
field description studied in [30] exactly coincides with the ABBM model. Long 
range elastic interactions lead to d = 3 being the upper critical dimension for 
earthquakes models [93,131,133,134], and playing the same role as the long 
range dipolar interactions in BN. 

It is interesting to mention that asymmetric avalanches are also observed in 
other crackling systems. The dynamics of the Bak-Tang-Wiesenfeld model of 
sand piles [135] is characterized by asymmetric avalanches [52]. Leftward asym- 
metric pulses also appear in the intermittent behavior of systems undergoing 
a shear deformation that proceeds in bursts of dislocations activity [136]. Re- 
cently Laurson and Alava [137] analyzed a discrete dislocation dynamics model 
and observed, in agreement with the experiments, a leftward asymmetry in the 
avalanche shape, corresponding to an asymmetry in the average creation and 
annihilation rates of dislocations during the avalanche. They also find that the 
asymmetry decreases with the signal threshold, although the threshold value 
do not have visible effects on the avalanche statistics. Finally, a strong skew- 
ness in the avalanche shape has recently been observed in the acustic emissions 
during martensitic transformations [138]. 



7 Conclusions 

In this review we report on the present state of understanding of the 
Barkhausen effect in soft ferromagnetic materials, focusing on those with long 
range dipolar interactions. In particular, we collect several results on the mean 
field ABBM model for BN. This model, originally proposed on phenomenolog- 
ical basis, corresponds to a mean field description for an elastic magnetic wall 
in a disordered ferromagnet under the effect of an external driving field, and 
is extremely successful in reproducing the complex phenomenology observed 
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in the experiments, at least as long as universal properties are involved. 

The ABBM model consists of an effective equation of motion for the velocity 
of the avalanche front. This equation maps onto the Langevin equation for a 
biased random walk in a logarithmic potential with appropriate boundary 
conditions. This process is amenable to analytical treatment with standard 
methods in the theory of stochastic processes. Therefore this mapping has the 
advantage to provide several exact results, as well as a direct interpretation of 
a number of properties observed in the experiments in terms of properties of 
the corresponding stochastic process. Specifically, through this mapping, the 
power laws in the distributions of sizes and durations of avalanches are exactly 
derived in terms of distributions of return times to the origin, the continuous 
dependence of the exponents in the power law distributions on the driving field 
rate can be calculated and turns out to be related to the marginality of the 
logarithmic perturbation to a free random walk, the existence of a threshold in 
the driving field rate has a clear interpretation (thanks to a further mapping) 
in terms of recurrence properties of a free random walk, the existence of a cut- 
off in the power law distributions is ascribed to the bias term, and therefore is 
due to the presence of the demagnetizing field, and finally the average shape 
of the Barkhausen pulse which is related to the excursion of the process, and 
the power spectra of the noise can be computed explicitly. 

Although deviations from the ABBM predictions are observed in some ex- 
periments, the picture that emerges from this theory is quite satisfactory. 
The most relevant experimental feature that is not reproduced by ABBM is 
the leftward asymmetric shape systematically observed in experiments. This 
asymmetry is due to a non-local damping by the eddy currents, whose effect 
is neglected in the original ABBM model formulation, as it is in most domain 
walls models. This effect can be taken into account by including in the equa- 
tion of motion an inertial term with a negative effective mass for the interface, 
which leads to the definition of a generalized ABBM model. The asymmetry 
of the pulses depend on the duration and can be used to extract important 
information on the characteristic time of the underlying dynamics. 

A basic point to keep in mind, is that the whole theory developed so far 
models one single domain wall, while many interacting domain walls are in 
general involved in the magnetization process. It would therefore be interesting 
to understand how the presence of many domains modifies the present theory. 

Although the ABBM theory was derived as a model for Barkhausen noise, 
it turns out to be a paradigmatic theory for avalanche dynamics, and can 
give suggestions on how to model other crackling systems. We discussed some 
analogies with sheared granular matter and seismic events. These similarities 
point to the existence of deep connections between the physics of avalanches 
in different systems. Also the presence of non-universal effects might be com- 
mon to other crackling system. Asymmetric avalanche shapes similar to those 
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observed in BN appear in other crackling systems: in martensitic transfor- 
mations, in plastic deformations, in the Bak-Tang-Wiesenfeld model of sand 
piles, and in seismic events. In the last case, the leftward asymmetric average 
moments may originate, from the presence of stress overshoots. 
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